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Nonlinear  approximation  from 
differentiable  piecewise  polynomials 

Oleg  Davydov  and  Pencho  Petrushev  * 


Abstract 

We  study  nonlinear  n-term  approximation  in  Lp(M2)  (0  <  p  <  oo)  from  hierarchical 
sequences  of  stable  local  bases  consisting  of  differentiable  (i.e. ,  Cr  with  r  >  1)  piece- 
wise  polynomials  (splines).  We  construct  such  sequences  of  bases  over  multilevel  nested 
triangulations  of  R2,  which  allow  arbitrarily  sharp  angles.  To  quantize  nonlinear  n- 
term  spline  approximation,  we  introduce  and  explore  a  collection  of  smoothness  spaces 
(B-spaces).  We  utilize  the  B-spaces  to  prove  companion  Jackson  and  Bernstein  esti¬ 
mates  and  then  characterize  the  rates  of  approximation  by  interpolation.  Even  when 
applied  on  uniform  triangulations  with  well-known  families  of  basis  functions  such  as 
box  splines,  these  results  give  a  more  complete  characterization  of  the  approximation 
rates  than  the  existing  ones  involving  Besov  spaces.  Our  results  can  easily  be  extended 
to  properly  defined  multilevel  triangulations  in  Md,  d  >  2. 


1  Introduction 

Nonlinear  approximation  of  functions  in  dimensions  d  >  1  is  a  challenging  area,  especially  if 
one  moves  away  from  tensor  product  type  approaches  in  order  to  more  adequately  approxi¬ 
mate  functions  with  singularities  along  curves  and  with  other  anisotropies.  One  of  the  most 
natural  tools  for  approximation  is  piecewise  polynomials  over  triangulations,  and  a  funda¬ 
mental  problem  is  to  characterize  the  rate  of  nonlinear  approximation  in  Lp  (0  <  p  <  oo)  in 
terms  of  properly  defined  global  smoothness  conditions.  This  problem  is  disheartening  if  one 
allows  the  nonlinear  approximation  to  be  from  any  piecewise  polynomial  over  an  arbitrary 
triangulation.  The  difficulty  stems  from  the  highly  nonlinear  nature  of  piecewise  polynomials 
in  dimensions  d  >  1.  For  instance,  if  sq  and  s2  are  two  piecewise  polynomials  over  n  triangles 
in  R2  each,  then  Si  +  s2  is  in  general  a  piecewise  polynomial  over  many  more  than  n  (even 
>  n2)  pieces.  Therefore,  the  well-known  recipe  of  proving  Jackson  and  Bernstein  estimates 
and  then  applying  interpolation  is  useless. 

The  problem  becomes  even  harder  when  differentiable  piecewise  polynomials  are  needed, 
which,  for  instance,  is  the  case  for  numerous  practical  applications  of  surface  modeling  and 
for  the  conforming  finite  element  methods  for  higher  order  PDEs.  Moreover,  there  is  an 
intrinsic  demand  for  differentiability  of  the  approximating  tools  from  the  point  of  view  of  the 
nonlinear  approximation  theory  itself.  Indeed,  this  property,  together  with  local  reproduction 
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of  higher  degree  polynomials,  is  crucial  for  the  ability  to  represent  higher  order  smoothness 
spaces,  such  as  classical  Sobolev  or  Besov  spaces  in  regular  settings  (see  Theorem  2.15). 
The  desirable  differentiability  of  the  approximating  piecewise  polynomials,  however,  leads  to 
additional  difficulties  because  of  the  complicated  structure  of  spaces  of  multivariate  splines. 
For  example,  the  dimension  is  not  known  and  stable  local  bases  are  impossible  in  general 
already  for  the  space  of  all  piecewise  polynomials  of  degree  <  k  and  smoothness  r  >  1  with 
respect  to  a  finite  triangulation  of  a  polygonal  domain  in  R2  if  k  <  3 r  +  2  [10]. 

A  reasonable  alternative  to  “spline  approximation  with  free  triangulations”  is  to  consider 
nonlinear  n-term  approximation  from  hierarchical  sequences  of  spline  bases  over  multilevel 
nested  triangulations  of  Rd.  (For  the  sake  of  simplicity,  we  shall  restrict  ourselves  in  this 
article  to  the  case  d  =  2.)  To  explain  this  concept  more  precisely,  consider  a  sequence 
°f  partitions  of  R2  into  triangles  with  disjoint  interiors  such  that  each  level  Tm.  is 
a  refinement  of  the  previous  one  Tm-i ■  Let  T  :=  Umez^m-  We  imP°se  certain  mild  (and 
natural)  conditions  on  the  triangulations  which  prevent  them  from  deterioration  but  still 
allow  the  triangles  to  change  in  size,  shape,  and  orientation  quickly  when  moving  around  at 
a  given  level  or  through  the  levels.  In  particular,  triangles  with  arbitrarily  sharp  angles  may 
occur  at  any  location.  We  denote  by  Sk'r(Tm )  the  set  of  all  r-times  differentiable  piecewise 
polynomials  with  respect  to  Tm  of  degree  <  k.  Given  a  ladder  of  spaces 

•••  C  S-x  C  <S0  C  Si  C  •••,  Sm  C  Sk'r(Tm),  (1-1) 

and  bases  of  Sm,  m  G  Z,  we  set  <E>  :=  $7-  :=  |JmeZ  Using  the  standard  wavelet 
terminology,  we  can  describe  such  a  nested  sequence  of  spaces  with  bases  as  “spline  mul¬ 
tiresolution”  (or  “multiresolution  analysis”). 

Consider  now  the  problem  for  nonlinear  (n-term)  approximation  from  the  set  En  of  all 
piecewise  polynomials  of  the  form  s  =  Xq-i  GAP  where  ip$  G  $  may  come  from  different  lev¬ 
els  and  locations.  Once  a  particular  multilevel  triangulation  has  been  selected,  the  variety  of 
piecewise  polynomial  approximations  significantly  reduces.  However,  a  great  deal  of  flexibil¬ 
ity  is  retained,  and  the  problem  remains  highly  nonlinear.  For  instance,  thin  and  elongated 
basis  functions  are  allowed.  On  the  other  hand,  the  advantages  of  multilevel  approximation 
methods  can  be  exploited  in  full. 

Our  program  consists  of  the  following  basic  steps: 

1.  We  construct  hierarchical  sequences  of  bases  ($m)mez  on  multilevel  triangulations  sat¬ 
isfying  certain  requirements  of  local  regularity  allowing  anisotropically  shaped  triangles. 

2.  To  quantify  the  approximation  process,  we  introduce  and  develop  a  family  (library)  of 
smoothness  spaces  5“  ($7-)  depending  on  $7-  and  as  a  consequence  on  the  triangulation 
T.  We  call  them  B-spaces  since  they  have  some  resemblance  to  Besov  spaces.  So,  the 
idea  is  to  measure  the  smoothness  of  the  functions  using  a  family  of  space  scales  B®  ($7-) 
(which  vary  with  $7-)  instead  of  a  single  scale  of  smoothness  spaces  like  the  scale  of 
Besov  spaces. 

3.  We  develop  a  coherent  theory  of  nonlinear  n-term  approximation  from  $7-  based  on 
the  idea  of  proving  Jackson  and  Bernstein  estimates  and  interpolation. 
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4.  We  utilize  this  theory  in  the  development  of  algorithms  for  nonlinear  piecewise  poly¬ 
nomial  (spline)  approximation  which  capture  the  rate  of  the  best  approximation. 

The  logic  of  the  resulting  approximation  scheme  is  the  following:  Suppose  {$r}r  is  a 
collection  of  multilevel  sequences  of  (spline)  bases  as  above. 

(i)  For  a  given  function  /,  find  the  “right”  triangulation  T  :=  7}  such  that  /  exhibits 
the  most  smoothness  (sparsity  of  its  representation)  when  measured  via  the  scale  Ba(Q T). 

(ii)  Find  an  optimal  or  near  optimal  representation  of  /  using  <Fr.  (Note  that  $7-  is 
redundant,  i.e.,  linearly  dependent.) 

(iii)  Using  this  representation,  run  an  algorithm  for  n-term  Lp-approximation  which 
achieves  the  rate  of  the  best  n-term  approximation. 

Naturally,  the  first  step  presents  the  most  challenging  problem  in  this  scheme.  We  do 
not  have  a  completely  satisfactory  algorithm  for  this  step.  (Note  that  this  problem  has  a 
complete  and  efficient  solution  in  the  simpler  case  of  nonlinear  approximation  from  piecewise 
polynomials  over  dyadic  partitions,  see  [54].)  As  it  will  be  shown  in  this  article,  the  other 
steps  are  now  well  understood  and  have  complete  solutions. 

The  above  program  has  been  suggested  in  [38]  and  implemented  in  [38,  39]  in  the  cases 
of  approximation  from  discontinuous  piecewise  polynomials  and  continuous  piecewise  linear 
functions  (r  =  —1,  k  >  1,  and  r  =  0,  A:  =  2,  where  r  =  —1  corresponds  to  the  discontinuous 
case).  The  simplest  example  of  a  hierarchical  family  of  continuous  basis  functions  is  the  set 
of  all  Courant  elements  generated  by  a  multilevel  nested  triangulation  T,  that  is,  the  set  of 
all  piecewise  linear  and  continuous  functions  $7-  =  {p>e}  supported  on  the  cells  {0}  (each  6 
is  the  union  of  all  triangles  of  a  particular  level  Tm  attached  to  a  vertex),  see  [38]. 

In  the  present  article,  we  develop  the  theory  of  nonlinear  n-term  approximation  for  basis 
families  consisting  of  differentiable  piecewise  polynomials  (r  >  1).  The  construction  of 
such  basis  functions  suitable  for  application  is  hampered  by  the  fact  that  both  the  classical 
differentiable  finite  elements  [14]  and  the  earlier  polynomial  spline  basis  constructions  on 
arbitrary  triangulations  [1,  8,  16,  17,  18,  35,  36,  44,  48,  57]  are  difficult  to  use  for  our 
purposes,  see  Remark  4.12  and  the  discussion  in  §5.3.  The  stable  local  spline  bases  of  [27] 
can  in  principle  be  used  in  two  variables.  However,  all  other  arguments  of  our  article  are 
basically  “dimension  independent”,  and  we  refrain  here  from  treating  the  case  d  >  2  only 
for  the  sake  of  simplicity  and  clarity.  Therefore,  we  build  upon  the  nodal  spline  bases  of 
[22],  which  is  the  only  known  approach  that  produces  stable  local  bases  for  nested  spline 
spaces  on  general  triangulations  in  all  dimensions.  However,  these  bases  are  stable  only  for 
triangulations  satisfying  (in  M2)  the  minimal  angle  condition.  We  extend  the  construction 
of  [22]  to  a  wider  class  of  strong  locally  regular  triangulations,  see  §2  for  a  definition.  Note 
that  the  new  basis  functions  are  invariant  under  affine  transforms  (see  Remark  4.9).  In  the 
case  r  =  0  our  construction  reduces  to  the  classical  continuous  Lagrange  finite  elements  and 
is  valid  for  any  locally  regular  triangulation,  see  Remark  4.13. 

A  focal  point  of  our  development  is  the  characterization  of  nonlinear  n-term  approxi¬ 
mation  from  families  of  differentiable  spline  basis  functions,  including  the  development  of 
B-spaces,  proof  of  Jackson  and  Bernstein  estimates,  and  characterization  of  the  approxima¬ 
tion  spaces  by  interpolation  (see  §2  §3).  In  [39],  there  are  three  algorithms  developed  for 
nonlinear  n-term  approximation  in  Lp  (0  <  p  <  00)  from  Courant  elements.  These  can  be 
immediately  implemented  for  n-term  approximation  from  differentiable  spline  bases,  and  it 
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can  be  shown  similarly  as  in  [39]  that  they  achieve  the  rate  of  the  best  approximation.  We 
do  not  pursue  this  goal  here. 

The  B-spaces  from  the  present  article  can  be  viewed  as  a  generalization  of  the  “ap¬ 
proximation  spaces”  from  §3.4  of  [51]  (see  also  the  references  therein).  More  precisely,  in 
the  specific  setting  of  “quasi-uniform  partitions”  and  the  basis  functions  used  in  [51],  our 
B-spaces  coincide  with  the  approximation  spaces  of  [51]. 

The  theory  of  nonlinear  n-terrn  approximation  from  box  splines  (on  uniform  triangula- 
tions)  has  been  developed  in  [29]  (p  <  oo)  and  [30]  (p  =  oo)  (for  nonlinear  spline  approxi¬ 
mation  in  dimension  d  =  1.  see  [53]).  In  these  articles,  direct,  inverse,  and  characterization 
theorems  have  been  proved  utilizing  certain  Besov  spaces.  Even  in  this  case,  our  results 
which  utilize  B-spaces  (in  place  of  Besov  spaces)  are  more  complete  since  they  character¬ 
ize  nonlinear  n-terrn  box  spline  approximation  for  all  rates  of  approximation  while  in  the 
above-mentioned  articles  the  rate  is  restricted  by  the  Besov  smoothness  of  the  box  splines. 

There  is  an  apparent  connection  between  our  developments  here  and  multilevel  finite 
element  methods  for  PDEs,  see,  e.g.,[51].  Therefore,  it  seems  an  interesting  task  to  de¬ 
velop  finite  element  algorithms  for  solving  PDEs  which  achieve  the  rate  of  the  best  n-term 
approximation  of  the  solution. 

The  outline  of  the  article  is  the  following.  In  §2,  we  introduce  and  develop  the  B- 
spaees  needed  for  the  characterization  of  nonlinear  approximation  for  any  family  of  basis 
functions  with  certain  properties.  In  §3,  we  develop  the  general  theory  of  nonlinear  n- 
term  approximation  from  piecewise  polynomials,  where  the  global  smoothness  of  functions 
is  measured  by  means  of  our  B-spaces.  In  §4,  we  construct  hierarchical  sequences  of  bases 
consisting  of  differentiable  piecewise  polynomials.  In  §5,  we  review  a  number  of  alternative 
constructions  fitting  into  our  scheme,  based  on  box  splines  and  some  other  spline  bases  on 
special  triangulations.  The  final  §6  is  an  appendix  containing  some  of  the  proofs. 

Throughout  the  article,  we  use  the  following  notation:  L]^C(R2)  :=  C(M2)  and  L^M2)  :  = 
C0(K2)  :=  {/  G  C(K2)  :  lim,  fix)  =  0},  if  :=  40C(M2),  0  <  q  <  oo,  C  :=  C(M2), 
II  '  lie  :=  II  '  I|l«(s:2),  0  <  q  <  oo;  11^  denotes  the  set  of  all  algebraic  polynomials  in  two 
variables  of  total  degree  <  k.  For  any  OcK2,  1q  denotes  the  characteristic  function  of  Q 
and  |D|  denotes  the  Lebesgue  measure  of  Q.  Positive  constants  are  denoted  by  c,  ci, . . .  (they 
may  vary  at  every  occurrence),  a  ~  [3  means  C\Ct  <  /3  <  c2ft,  and  a  :=  f3  or  [3  =:  a  stands 
for  “ft  is  by  definition  equal  to  /?”. 


2  B-spaces  generated  by  spline  multiresolution 

In  the  present  section,  we  introduce  and  explore  the  smoothness  spaces  we  need  for  the  char¬ 
acterization  of  nonlinear  n-term  spline  approximation  generated  by  families  of  differentiable 
basis  functions  over  multilevel  nested  triangulations. 

2.1  Triangulations 

In  our  development,  we  utilize  three  types  of  multilevel  nested  triangulations.  We  shall 
call  each  of  them  simply  a  triangulation ,  although  such  a  triangulation  does  not  form  a 
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single  partition  of  R2  but  rather  an  infinite  nested  family  of  partitions  (each  of  them  is  a 
triangulation  of  R2  in  the  more  commonly  used  sense). 

Let  T  =  Ume^Tm  be  a  set  of  closed  triangles  in  R2  with  levels  Tm ,  m  G  Z.  Denote  by 
Vm  the  set  of  all  vertices  (nodal  points)  of  triangles  from  Tm  and  set  V  :=  |JmeZ  Vm.  We  say 
that  T  is  a  triangulation  of  R2  if  the  following  conditions  are  fulfilled: 

(a)  Every  level  Tm  is  a  set  of  triangles  with  disjoint  interiors  which  cover  R2:  R2  = 

U 

(b)  The  levels  {Tm)mei  °f  T  are  nested ,  i.e. ,  Tm+i  is  a  refinement  of  Tm  obtained  by  splitting- 
each  A  G  Tm  into  subtriangles  with  disjoint  interiors  called  children  of  A. 

(c)  Each  triangle  A  G  Tm  has  at  least  two  and  at  most  M0  children  in  Tm.+\,  where  M0  >  2 
is  a  constant  independent  of  m. 

(d)  No  hanging  vertices  condition:  No  vertex  of  any  triangle  A  G  Tm  lies  in  the  interior  of 
an  edge  of  another  triangle  from  Tm- 

(e)  The  valence  Nv  of  each  vertex  v  G  Vm  (the  number  of  triangles  A  G  Tm  which  share  v 
as  a  vertex)  is  <  N0,  where  N0  is  a  constant. 

(f)  For  any  compact  K  C  R2  and  any  fixed  m  G  Z,  there  is  a  finite  collection  of  triangles 
from  Tm  which  cover  K. 

Note  that  any  two  triangles  in  T  either  have  disjoint  interiors  or  one  of  them  contains  the 
other.  In  particular,  A'  G  Tm+ 1  is  a  child  of  A  G  Tm  ( rn  G  Z)  if  and  only  if  A'  C  A.  If  A 
and  A'  are  two  different  triangles  in  T  and  A'  C  A,  then  we  say  that  A  is  an  ancestor  of 
A',  while  A'  is  a  descendant  of  A. 

Locally  regular  triangulations.  We  call  a  triangulation  T  =  Umez  7m  a  locally  regular 
triangulation  of  R2,  or  briefly  an  LR-triangulation ,  if  T  satisfies  the  following  additional 
conditions: 

(g)  There  exists  a  constant  1/2  <  p  <  1  such  that  for  each  A  G  T  and  any  child  A'gT 
of  A, 

(1  P)  I A  |  <  |  A'|  <  p|  A|.  (2.1) 

(hi)  There  exists  a  constant  0  <  <  1  independent  of  m  such  that  for  any  A',  A"  G  Tm 

(m  G  Z)  with  a  common  edge, 

A  <  | A'|/| A"|  <  S^1.  (2.2) 

By  (e),  it  follows  that  for  any  A',  A"  G  Tm  with  at  least  one  common  vertex,  (2.2)  holds 
with  <5i  replaced  by  h|Vo'/2. 

Strong  locally  regular  triangulations.  We  call  a  triangulation  T  =  Um&  7m  a  strong 
locally  regular  triangulation  of  R2  or  briefly  an  S LR-triangulation ,  if  T  satisfies  (2.1)  and  the 
following  condition  that  replaces  (2.2): 
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(h2)  There  exists  a  constant  0  <  82  <  1/2  such  that  for  any  A',  A"  G  Tm  ( m  G  Z)  sharing 
an  edge, 

|conv(A,U  A")|/| A'|  <  S^1,  (2.3) 

where  conv  (G)  denotes  the  convex  hull  of  G  C  M2. 

Obviously,  (2.3)  implies  (2.2)  with  A  =  82.  Therefore,  each  SLR-triangulation  is  an  LR- 
triangulation. 

Regular  triangulations.  By  definition,  a  triangulation  T  =  \Jrnev  Tm  is  called  a  regular 
triangulation  if  T  satisfies  the  following  condition: 

(h3)  There  exists  a  constant  f3  =  /3(T)  >  0  such  that  the  minimal  angle  of  each  triangle 
A  G  T  is  >  f3. 

Next,  we  make  a  few  remarks  which  will  help  understand  better  the  nature  of  the  trian¬ 
gulations  that  we  utilize. 

(i)  For  each  of  the  three  types  of  triangulations  there  is  a  number  of  constants  that 
are  assumed  fixed.  In  what  follows  we  refer  to  them  as  parameters.  Thus  the  parameters 
of  an  SLR-triangulation  are  M0,  N0,  p,  and  82.  Notice  that  because  of  (2.1),  we  can  set 
M0  :=  1/(1  —  p)  and  remove  M0  from  the  list  of  parameters.  However,  this  would  tend  to 
obscure  the  actual  role  of  p  and  M0. 

(ii)  It  is  a  key  observation  that  the  collection  of  all  SLR-triangulations  with  given  (fixed) 
parameters  is  invariant  under  affine  transforms.  The  same  is  true  for  LR-triangulations. 

(iii)  It  is  easy  to  see  that  (2.3)  is  equivalent  to  the  following  condition  introduced  in  [38]: 
Affine  transform  angle  condition:  There  exists  a  constant  fi  =  /3(T),  0  <  :i  <  7t/3,  such 

that  if  A0  G  Tm,  m  G  Z,  and  A  :  R2  — >■  M2  is  an  affine  transform  that  maps  A0  one-to-one 
onto  an  equilateral  reference  triangle,  then  for  every  A  G  Tm  which  has  at  least  one  common 
vertex  with  A0,  we  have 

minangle  (A(A))  >  /3,  (2.4) 

where  A(A)  is  the  image  of  A  by  the  affine  transform  A. 

The  equivalence  of  the  two  conditions  follows  easily  from  the  obvious  but  important  fact 
that  both  conditions  are  invariant  under  affine  transforms. 

Note  that  we  prefer  to  use  (2.3)  rather  than  (2.4)  in  the  definition  of  SLR-triangulations 
in  this  article  since  the  constant  82  appears  naturally  when  estimating  norms  of  the  basis 
functions  constructed  in  §4  (see  (4.8))  and  also  (2.3)  is  easier  to  verify  in  practical  situations. 

(iv)  As  we  have  already  mentioned,  every  SLR-triangulation  is  an  LR-triangulation  but 
the  converse  statement  is  not  true.  Also,  every  regular  triangulation  is  an  SLR-triangulation 
but  not  the  other  way  around.  Counterexamples  are  given  in  [38]. 

(v)  The  maximal  angle  (MA)  condition 

7 r  —  max  angle  (A)  >  /3  >  0,  A  G  T,  (2.5) 

known  from  the  finite  element  method  [2]  is  totally  different  from  our  conditions  of  regular¬ 
ity.  It  is  easy  to  see  that  there  are  SLR-triangulations  that  do  not  satisfy  MA,  and  there 
are  triangulations  that  satisfy  MA  and  fail  to  be  locally  regular.  As  we  shall  see  below 
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(Example  4.7),  our  construction  of  stable  differentiable  basis  functions  does  not  extend  to 
triangulations  satisfying  the  MA  condition  but  failing  to  be  SLR. 

(vi)  The  rate  of  change  of  the  size  of  the  elements  (|  A|,  min  angle  (A),  and  diarn(A))  of  a 
triangle  A  G  T  as  A  moves  away  from  a  fixed  triangle  A*  G  T  for  different  types  of  triangu¬ 
lations  T  is  explored  in  [38] .  We  shall  briefly  discuss  this  issue  for  SLR-triangulations  which 
are  the  most  important  type  of  triangulations  for  the  present  article.  An  SLR-triangulation 
T  may  have  an  equilateral  (or  close  to  such)  triangle  A*  at  any  level  Tm  with  descendants 
Ai  D  A2  D  •  ■  •  such  that  min  angle  (Aj)  — >■  0  as  j  — >■  oo,  and  also  a  sequence  (A')°L0  C  Tm 
with  A'0  =  A  and  A  '•  n  A'+1  ^  0  (j  =  0, 1, . . .)  such  that  min  angle  (A')  — >■  0.  Conditions 
(2.1)  and  (2.3)  suggest  geometric  rates  of  change  of  |A|,  min  angle  (A),  and  diarn(A)  as 
A  G  Tm,  moves  away  from  a  fixed  A*  G  Tm,-  In  fact,  the  rate  of  change  is  a  power  of  the 
minimal  number  of  edges  connecting  A  and  A°,  see  [38]. 

(vii)  We  shall  need  to  know  what  happens  with  the  levels  Tm  of  a  triangulation  T  as 

m  — >■  — oo.  By  Lemma  2.1  from  [38],  for  each  LR-triangulation  T  there  exists  a  finite  cover 
T-oo  °f  R2  such  that  either  T-oo  =  {R2}  or  T-oo  =  (A'x  )J\.  <  N0,  where  each  A^  is 

an  infinite  triangle ,  i.e.,  the  set  of  all  points  on  and  between  two  rays  which  are  not  collinear 
and  have  a  common  beginning.  Moreover,  in  the  second  case,  the  infinite  triangles  ( Afi  N, 
have  a  single  common  vertex  and  disjoint  interiors,  and  also  each  triangle  A  G  T  and  all  its 
ancestors  are  contained  in  an  infinite  triangle  A^  G  T-oo- 

For  more  details  about  multilevel  triangulations,  see  [38]. 

Some  additional  notation  and  preliminaries.  We  denote  by  [tq,  vfi  the  interval  (straight 
line  segment)  with  endpoints  v \ ,  v->  and  by  |e|  the  length  of  e  =  [tq,t>2].  Furthermore,  we 
let  [vi,V2,vz\  denote  the  triangle  with  vertices  tq,t>2,^3,  and  let  |A|  denote  the  area  of 
A  =  [tq,  v*h  tq].  Throughout  the  article,  we  assume  that  the  vertices  tq,  rq,  U3  of  any  triangle 
[uj,  u2,u3]  are  ordered  counterclockwise. 

For  a  triangle  A  G  Tm  ( m  G  Z),  we  define  level(A)  :=  m. 

For  any  vertex  v  G  Vm,  we  let  star  (v)  =  star1(u)  denote  the  star  of  v,  i.e.,  the  union  of 
all  triangles  A  G  Tm  attached  to  v.  Moreover,  for  each  t>  2,  we  denote  by  star  fiv)  the 
union  of  starl_1(u)  and  the  stars  of  the  vertices  of  staA_1(u).  (Note  that  star  fiv)  depends 
also  on  the  level  m.  but  we  do  not  indicate  this  in  the  notation  since  it  is  always  clear  from 
the  context  what  level  is  meant.)  We  also  set 

:=  U{starf(u)  :  v  G  Vm,  A  C  starve)},  A  G  Tm-  (2.6) 

It  is  easy  to  check  that  Qfi  :=  U{star  2£_1(u)  :  v  is  a  vertex  of  A},  A  G  Tm- 

It  is  readily  seen  that  there  exists  a  constant  c*  =  <f{N0 ,  l)  <  Nq  such  that 

#{A  G  Tm  :  A  C  star£(u)}  <  c*,  v  G  VTO,  (2.7) 

and  hence  there  exists  a  constant  c**  =  c**(N0,£)  <  3 c*(N0,  2i  —  1)  —  5  <  3Aq£_1  such  that 

#{AG7]„:Acfli,}<A,  A'  g  Tm-  (2.8) 

We  denote  by  £m  the  set  of  all  edges  of  triangles  of  Tm  and  set  £  :=  U mez^m-  We  let 
star(e)  denote  the  union  of  the  two  triangles  attached  to  e  G  £m. 
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(2.9) 


For  future  use,  we  state  the  following  inequality: 

OO 

(lA,|/lAD7  <  ^P’1  =  C(P.7)  <  °o,  A'  e  T,  7  >  0, 

Aer,  Ad  A'  i=o 

which  is  immediate  from  the  properties  of  LR-triangulations  (|A'|  <  p|A|  if  A'  is  a  child  of 
A). 

2.2  Basis  functions:  The  general  setting 

Let  T  =  Umez  Tn  he  a  locally  regular  (or  better)  triangulation.  For  m  E  Z,  r  >  0,  and 
k  >  1,  we  denote  by  =  Sk'r(Tm)  the  set  of  all  r  times  differentiable  piecewise  polynomial 
functions  of  degree  <  k  over  Tm ,  be.,  s  E  if  and  only  if  s  E  Cr  (R2 )  and  s  =  ^a'-Fa 

with  Pa  €  LR.  Naturally,  will  denote  the  set  of  all  piecewise  polynomials  of  degree 

<  k  over  Tm  which  are,  in  general,  discontinuous  across  the  edges  from  £m. 

We  assume  that  for  each  m  E  Z  there  is  a  subspace  Sm  of  Skf  ( r  >  0,  k  >  2)  and  a 
family  <Fm  =  {<pg  :  9  E  ©m}  C  Sm  of  basis  functions  satisfying  the  following  conditions: 

1.  n^.  C  Sm  for  some  1  <  k  <  k  (k  independent  of  m). 

2.  Sm  C  Sm. |_i  (m  E  Z). 

3.  For  any  s  E  Sm  there  exists  a  unique  sequence  of  real  coefficients  a(s)  =  {cie(s)) oeem 
such  that 

s  =  ^2  ao(s)£o- 
oeem 

(Thus,  is  a  basis  for  Sm  and  {ae{-))ee&m  are  the  dual  functionals.) 

4.  For  each  9  E  Qm  there  is  a  vertex  v  =  vg  E  Vm  such  that 

supp  lPq  C  staA(u)  =:  Eg, 

lbelUoo(K2)  =  WveWLooiEg)  <  Mi, 

|«o('S)|  <  M2||s||Loo(i^;),  s  E  Sm, 

where  t  >  1  and  M1,  M2  are  positive  constants,  all  independent  of  9  and  m. 

Let 

$  :=  U  Tm  and  0  :=  |J  0TO. 

We  shall  refer  to  r,  k,  k ,  t,  M} ,  and  M2  as  parameters  of  <F. 

A  simple  example  of  a  family  of  basis  functions  satisfying  the  above  conditions  is  the  set 
of  well-known  Courant  elements  (continuous  piecewise  linear  basis  functions,  r  =  0,  k  =  2) 
associated  with  T  (see  [38]).  Concrete  constructions  of  differentiable  basis  functions  (r  >  1) 
will  be  discussed  below  in  Sections  4  5. 

Although  0  and  0TO  (m  E  Z)  are  simply  index  sets,  in  the  case  of  Courant  elements, 
0  can  be  identified  as  the  set  of  all  cells  (supports  of  basis  functions).  As  we  shall  see  in 
Sections  4  5,  in  general,  several  basis  functions  of  <Fm  may  have  the  same  support.  However, 
the  supports  of  only  <  constant,  of  them  may  overlap: 


(2.10) 

(2.11) 

(2.12) 
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Lemma  2.1.  There  is  a  constant  L  depending  only  on  k,  (',  and  N0  such  that  for  any  A  G  Tm 
(m  G  Z), 

if  {9  e  em 

where  Eg  is  defined  in  (2.10). 

Proof.  We  have  by  (2.10)  and  (2.8) 

if  {9  G  0m  :  A  C  Eg}  < 


< 

We  shall  frequently  use  the  equivalence  of  different  norms  of  polynomials  as  stated  in  the 
following  lemma  (see  also  [38]). 

Lemma  2.2.  Let  P  G  1U,  A:  >  1,  and  0  <  p,  q  <  oo. 

(a)  For  any  triangle  A  C  M2,  ||P||l3J(a)  ~  |  A|1/p_1/|?||P||i<;(A)  with  constants  of  equiva¬ 
lence  depending  only  on  p,  q,  and  k. 

(b)  If  A  and  A'  are  two  triangles  such  that  A'c  A  and  |A|  <  ci|A'|,  then  ||P||lp(A)  < 
H  P\  with  c  =  c{p,  A:,  cfi). 

(c)  If  A'  and  A  are  two  triangles  such  that  A'c  A  and  |A'|  <  C2 1 A  |  with  0  <  C2  <  l, 
then  ||P||L  (a)  <  c||P||Lp(A\A')  ~  |^|1//p_1//,J||-P|U<;(A\A')  with  constants  depending  oiily  on  p, 
q,  k,  and  c2. 

By  (2.2)  and  (2.7),  \Eg\  &  |A|  if  A  C  Eg ,  A  G  Tm,  and  9  G  0m.  Using  this  and 
Lemma  2.2,  we  obtain  that,  for  0  <  p,  q  <  oo, 

lls'lliP(£;#)  ~  \Eo\^p  1//,5,|ls'll Lq(Ee),  s  G  Sm,  9  G  0m,  (2.14) 

where  the  constants  of  equivalence  depend  on  p,  q.  k,  and  hi.  In  particular,  we  shall  need 
(2.14)  with  s  =  pg,  when  it  takes  the  form  \\pg\\p  ~  lEgf/P-^Wpg  ||g,  in  view  of  (2.10). 

Lemma  2.3.  The  bases  <bm  are  Z^-stable  for  all  0  <  q  <  oo.  That  is,  if  g  :=  hope, 

where  {bg)geem  is  an  arbitrary  sequence  of  real  numbers,  then 

mi,  » ( E  iiw»ii;)1/s- 

oeGm 


:  Eg  D  A}  <  L,  (2.13) 


<  dim Sk'  l{Tm)\^ 

(y  ^  ^^{A7  e  Tm  ■  A'  C  TllA} 


Moreover,  for  any  7  G  M  and  0  <  r  <  oo, 

(  E  (iApiisik<Ai)T)1/T « ( E  (|ai7ii»«p«ii,)t)1/t  (2.i5) 

A  ETm  OeGm 

where  the  constants  of  equivalence  are  independent  of  rn  and  g.  In  the  case  q  =  00  (or 
t  =  00)  the  Iq-norm  (. iT-norm )  above  is  replaced  by  the  sup-norm  as  usual. 
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Proof.  We  have  to  prove  only  (2.15),  since  the  first  statement  of  the  lemma  then  follows 
with  7  =  0  and  r  =  q.  For  each  A  G  Tm,  we  have  by  (2.10), 

I '/I  /  =  II  boVe\\q  <c  \\boVo\\q- 

6e&m,E0DA  oeem,EeD  a 

Therefore,  by  Lemma  2.1  and  (2.7), 

E  (iAriisik,(A>)T  <  A  E  <i£«riM,>T 

A eTm  AeTm  Oe&m  ,  EgDA 

<  o  E  (i£»riiv»iid. 

oe&m 

In  the  other  direction,  since  $  is  a  basis  of  Sm  and  g  G  Sm .  we  have  bo  =  ag(g ),  9  G  0m, 
and  hence,  by  (2.12),  (2.14),  and  (2.11), 

\\be^o\\q  =  \\^e{g)^o\\q  <  <]  g\  I  N;ft»)|  rll\  q  <  <]  <l\  I  /•«  '  '' 

<  c\\g\\LqiEe)  <c  \u\  / 

&£Tm,AcEe 

Since  \Eg\  ~  |A|  if  A  G  Tm  and  A  C  Eg ,  we  have,  by  (2.7)  and  Lemma  2.1, 

E(Wii%Ar  <  A  E  dAriisik(A)r 

Oe&m  Oe&m  AeTm,  ACBfl 

<  c  E  (lAPllslli,W))T-  ■ 

A  eTm 

Local  polynomial  approximation  is  an  important  tool  in  spline  approximation.  For  a 
function  /  G  Lq(G ),  Gel2,  we  denote  by  Ek(f ,  G)q  the  error  of  the  best  Z^-approximation 
to  /  on  G  from  11^  and  by  04(/,  G)q  the  Adli  local  modulus  of  smoothness  of  /  on  G : 

Ek{f,G)q  :=  inf  \\f  -  P\\Lq(G),  Uk(f,  G)q  :=  sup  \\Akh(f, -)||i,(G)- 
PeUk  hew 

Whitney’s  theorem  gives  an  important  relation  between  these  two  quantities:  If  /  G  Lq(G ), 
0  <  q  <  oo,  where  G  =  A  is  an  arbitrary  triangle  or  G  =  a  with  A  G  T,  T  an  SLR- 
triangulation,  then 

Ek(f.G)q<c^.(f.G)q,  (2.16) 

where  c  =  c(q ,  k)  if  G  =  A  and  c  =  c(q,  A:,  62)  if  G  =  12a  (<E  is  from  (2.3)).  For  a  proof  of 
this  estimate,  see,  e.g.,  the  appendix  of  [38].  Note  that  this  estimate  holds  for  much  more 
general  regions  G,  but  then  the  constant  c  =  c(G)  may  become  hard  to  control. 

For  0  <  q  <  00  and  a  triangle  A,  we  let  Pa,9  :  Lq(A)  — >■  II/.  be  a  projector  such  that 

11/  -  P%(/)||l,(a)  <  cEk(f ,  A)q  for  /  G  Lq( A).  (2.17) 

Note  that  Pa,9  can  be  realized  as  a  linear  projector  if  q  >  1.  For  instance,  one  can  utilize 
the  averaged  Taylor  polynomial.  Namely,  suppose  A0  is  an  equilateral  reference  triangle 
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and  A  is  an  affine  transform  mapping  A  onto  A0.  Let  now  P(g)  G  11^  be  the  averaged 
Taylor  polynomial  of  the  function  g  :=  f  o  A-1  (the  composition  of  /  with  A-1)  over  the 
disc  B  inscribed  in  A0  (see,  e.g.,  §4A  of  [12]).  Clearly,  P  :  Lq(B)  — >■  11/,  is  a  linear  operator, 
\\P(g)\\Lq(B)  <  c\\g\\Lq(B)  (Q  >  1),  and  P  is  a  projector,  i.e.,  P(Q)  =  Q  for  Q  G  Uk.  From 
these  properties  of  P,  it  follows  that  for  an  arbitrary  Q  G  11*., 

h  -  P(g)\\Lq(A0)  <  \\g  -  Q\\Lq(/\0)  +  \\Q  -  P{g)\\Lq(A0) 

<  h  -  Qlk(Ao)  +  c\\P(g  -  Q)\\Lq(B)  <  c\\g  -  Q[|l,(a0) 

which  implies  \\g  —  P(g)\\Lq(A0)  A  cEk(g,  A0)?.  Substituting  back,  one  easily  obtains  \\f  — 
(P  o  .  1  )(  /')|  /  A  cEk(f,  A)q.  Finally,  we  set  Paj9  :=  P  o  A,  which  is  the  desired  linear 
projector  of  Lq(A)  into  II/.  . 

Note  that  Pa,§  cannot  be  realized  as  a  linear  operator  if  0  <  q  <  1  (otherwise,  we  would 
be  able  to  construct  a  nonzero  bounded  linear  functional  on  Lq). 

We  define  a  linear  operator  Qm  :  Sk~1(Tm)  — >■  Sm  as  follows.  For  each  6  G  ©m,  let 
\e-Sk’-\Tm)\  Ee  —>  R  be  a  linear  functional  such  that 

s  ^  and 

\Mf)\<M2\\f\\LooiEeh  f  e  Sk'-\Tm)W 

Such  linear  functional  always  exists  by  the  Hahn-Banach  theorem.  We  set 

Qm(s)  :=  ^2  'WSI-E«)^A  S  G  sk'  1  ( Irn )  •  (2.18) 

oeem 

Clearly,  Qm(s )  =  s  if  s  G  S„ and  thus  Qm  is  a  linear  projector  of  Sk~l(Tm )  into  <Sm . 
Lemma  2.4.  For  any  s  G  Sk~1(Tm),  0  <  q  <  oo,  and  A  G  Tm, 

|  Q m ( -s' ) |  i.9{/\ ;  A  cllsllL,J(n^)i  (2.19) 

with  a  constant  c  independent  of  m,  A,  and  s. 

Proof.  By  Lemma  2.2  and  (2.14),  we  have 

IMma)  <  ci|A|1/<?||v90||ioo(A)  <  CXM il A|1/<?, 

INUcAPO  <  c2\A\~L/q\\s\\Lq(Ee), 

where  C\  and  c2  depend  only  on  q  and  k.  Therefore, 

||<2m(s)||&i(A)  =  y  <  C  |Ag(g|^)|  |  r~n\  /„;a; 

ee®m  q(  '  ee@m 

Ac  Eq  Ac  Eg 

<  c  5]  llsllioc,(^.)|A|1/9  <  C  ||s||l9(S#)  <  c||s||£fl(n^. 

0  G  ©m  #G  ©m 

AC  Eg  ACE0 


11 


We  now  extend  Qrn  to  Ll°c(R2),  0  <  q  <  oo.  Let  P^,q  :  Lq( A)  — >■  n*.  be  a  projector 
satisfying  (2.17)  (linear  if  q  >  1).  We  define 

PmM )  ■■=  E  1A  •  PA1<?(/)  for  /  G  (2.20) 

AeTm 

which  is  a  projector  of  L^oc  into 
We  put 

<9rn,,(/)  :=  Qm(PmM))  for  /  G  Lf,  (2.21) 

which  is  evidently  a  projector  of  L!°c  into  <Sm  (linear  if  q  >  1  and  all  P^,q  are  linear). 

We  next  show  that  Qm,q  provides  a  good  local  Z^-approximation  from  Sm.  We  let  Sa (f)q 
denote  the  error  of  Lq ( Qfi ) -approxi rnati on  from  Sm .  i.e. , 

Sa (f)q  :=  inf  \\f  -  s||M^),  A  G  Tm.  (2.22) 

Thus,  Sa (f)q  is  the  error  of  approximation  to  /  from  restrictions  to  0/  of  functions  from 
Sm .  which  is  not  necessarily  the  same  as  the  approximation  by  all  r  times  differentiable 
piecewise  polynomials  of  degree  <  k  defined  only  on  even  if  Sm  coincides  with  Sffi . 
However,  since  C  Sm .  Sa (f)q  does  not  exceed  the  error  of  Lq ( ) -approxi mat i on  to  / 
from  polynomials  of  degree  <  k. 

Lemma  2.5.  If  f  G  Ll°c(R2),  0  <  q  <  oo  (/  G  C  if  q  =  oo),  then 

11/  -  Qm,,(f)\  /  <  cSA(/)9,  AG/m  {m  G  Z), 

with  c  independent  of  f,  m,  and  A. 

Proof.  Let  s. &  G  Sm  be  such  that  \\f  —  sa|Jl,(o^)  —  c§a (f)q-  Using  the  properties  of  Qm 
(see  Lemma  2.4),  we  find 

|  .1  Qm.<i[f)\  /  '  :  =  \\f  ~  Qm{Pm,q{f))  ||%(A| 

—  c||  /  •§  A  I  |i(J(  A)  A  c|  -s' a  Qm{Pm,q{f))  ||l9(A) 

A  cSa {f)q  A  c||Qm(SA  Pm.,q{f) )  ||l9(A) 

<  cSa(/)„  I  <]  •‘-•A.  p,,,„(f)\  / 

<  cSA(/),  a  c||/  -  sa||l,(^a)  a  c||/  -  Pm.„(./')|  /  ,((o<a; 

<  cSA(/)9.  ■ 

Lemma  2.6.  (a)  If  f  G  L^0C(R2)?  0  <  q  <  oo,  then  for  every  compact  K  C  R2, 

[|/ _  Qm,g(/)||Lapr)  -A  0  as  m  — >  oo.  (2.23) 

(b)  If  f  G  L?(M2),  0  <  <7  <  oo,  then 

\\f  -  Qm,q{f)\\Lq(wt?)  -A  0  as  m  — >■  oo.  (2.24) 

For  the  proof  of  this  lemma,  we  need  the  following  result: 
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Lemma  2.7.  If  T  is  an  LR-triangulation,  then  for  each  triangle  A*  G  T 

max{diam(A)  :  A  G  Tm ,  A  C  A°}  —y  0  as  m  — >■  oo.  (2.25) 

Proof.  Let  mo  :=  level(A°).  We  set  dm  :=  max{diam(A)  :  A  G  7m,  A  C  A°}.  Since 
( dm)m=m0  is  non-increasing,  it  suffices  to  show  the  existence  of  a  subsequence  tending  to 
zero.  Let  e  be  an  edge  of  a  triangle  A  G  Tm,  A  C  A°.  If  it  is  also  an  edge  of  a  child  of  A, 
then  the  valence  of  at  least  one  of  the  two  endpoints  of  e  will  increase  by  one  at  level  m  +  1. 
(Recall  that  there  are  always  at  least  two  children,  so  that  a  child  and  a  parent  cannot  be  the 
same  triangle.)  Therefore,  e  will  be  subdivided  at  least  once  after  at  most  S  :=  2 (N0  —  3)  + 1 
steps  of  refinement.  By  (2.1),  it  readily  follows  that  any  edge  e'  obtained  by  subdividing  e 
satisfies  |e'|  <  p\e\  <  pclm. 

We  call  an  edge  of  a  descendant  of  A®  a  cutting  edge  for  A*  if  one  of  its  endpoints 
is  a  vertex  of  A*  and  the  other  lies  in  the  interior  of  the  opposite  edge  of  A*.  Since  all 
cutting  edges  must  emanate  from  the  same  vertex  of  A*,  there  are  totally  no  more  than 
M  :=  N0  —  3  such  edges  for  A°.  Therefore,  no  new  cutting  edges  for  A*  will  be  created  at 
levels  m  >  m-o  +  M.  (It  is  easy  to  see  that  as  soon  as  no  new  cutting  edges  are  created  at  a 
level  m,  they  cannot  be  created  on  any  further  level.)  Using  this  and  the  above  observation, 
we  conclude  that  there  will  be  no  cutting  edges  at  levels  m  >  m0  +  M  +  S  since  they  all 
will  be  subdivided.  Therefore,  each  edge  e  inside  A°  at  these  levels  is  either  a  proper  part 
of  an  edge  of  A4  or  has  both  of  its  endpoints  in  the  interiors  of  two  different  edges  of  A*,  or 
it  has  at  least  one  endpoint  in  the  interior  of  A*.  In  all  cases,  condition  (2.1)  ensures  that 
\e\  <  M„o  •  which  implies  dmi  <  pdmo ,  where  m i  =  m0  +  M  +  S  +  1.  It  is  clear  now  that 
there  is  an  increasing  sequence  {m*.  }^=l  such  that 

dm j.  U  P  dmg  y  0  as  k  y  oo, 
which  completes  the  proof.  ■ 

Proof  of  Lemma  2.6.  (a)  By  condition  (f)  on  triangulations,  it  suffices  to  prove  the  lemma 
for  K  =  A*,  an  arbitrary  triangle  from  T-  Bv  Lemma  2.7, 

niax{diani(O^J  :  A  G  Tm,  A  C  ^0  as  m  ~ ^  °°-  (2.26) 

Case  1:  q  <  oo.  Fix  £  >  0.  In  view  of  (2.26),  there  exists  a  piecewise  constant  function 
S£  of  the  form 

S£  =  ^  ca1a,  h>.  >  level  (A4), 

AeTme, 

such  that 

\\f  ~  S£\\LqiniAo)  <  £  (2-27) 

(choose  first  g  G  C(QfA:/)  so  that  \\f  —  g\\Lq(nl  0)  <  -/2  and  then  choose  S£  so  that  \\g  — 

The"  Qm.«(Se)  =  QAS ,)■ 

We  have,  for  m  >  w . . 

11/  “  Qm.,/U  )\  /.,,(/  i  <  C\\f  ~  S.  |  +  c\\S£  —  Qm,q{S£)\\Lq(Ao) 

+  c\\Qm(S£-PmM))\\Lq(Ao).  (2.28) 
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For  the  third  term  above,  we  have 

\\Qm{Ss  —  Pm,q{f))  ||i9(A«)  <  c\\^e  ~  Ptn.qif  )\\ Lq(QfA<>) 

<  c\\f  ~  Se  Hl9(Q^)  +c\\f  -  Pm,q{f)  11^(0^) 

<  4f  -  S*\\Lq(^.)  <  C£,  (2.29) 

where  we  used  Lemma  2.4  and  that  \\f  -  Pm,q(f)\\Lq(rtA<>)  <  c|| /  -  (m  >  me),  by 

(2.17). 

It  remains  to  show  that  ,(.S',)|  r.,'A<>)  <  ce  for  sufficiently  large  m.  Denote 

by  G  the  union  of  the  edges  of  all  triangles  A  G  Tms  such  that  A  C  A°,  and  by  G$  :  = 
{rc  G  R2  :  dist(x,G)  <  5}  the  ^-neighborhood  of  G.  Clearly,  there  exists  S  >  0  such  that 

\\Se\\Lq(Gs)  <  £- 

By  (2.26),  there  exists  m i  >  rn~  such  that  diam(QA)  <  <5  for  all  triangles  A  G  Tm 
(m  >  mi)  such  that  A  C  A*  and  QA  fl  G  /  0.  Since  II j  C  <STO,  Q„,(S. )  ,•  =  Ss |a  if 
S£\ =  constant..  Using  this,  we  obtain  by  Lemma  2.5 

114  -  Q„l,,(Se)||lt(A.)  <  c(  5]  SA(S«)j) 1/5  <  clKII^  <  a. 

A  eTm,fFAnG^0 

We  substitute  this  estimate  together  with  (2.27)  and  (2.29)  in  (2.28)  to  obtain 

\\f-QmM)\\L  q(A<>)  <  ce  for  m  >  m.\. 

This  implies  (2.23)  if  q  <  oo. 

Case  2:  q  =  oo.  We  have,  by  Lemma  2.5  and  the  fact  that  IIi  C  Sm . 

||/  -  Qm,q{f)  ||ioo(A°)  <  ||/  -  C\ 

Now  the  result  follows,  using  (2.26)  and  the  fact  that  /  is  uniformly  continuous  on 
Part  (b)  of  the  lemma  is  immediate  from  part  (a).  ■ 

We  denote  <S_X  :=  f)meZSm.  As  we  already  mentioned,  there  are  only  two  possibilities 
for  71  oo:  71  oo  =  {R2}  or  71  oo  =  ( A/  A  .  AA  <  A0,  where  { A/^ }  are  infinite  triangles  with 
disjoint  interiors  and  a  common  vertex  which  cover  R2.  If  T-oo  =  {R2},  then  obviously  R2 
is  the  union  of  a  sequence  of  nested  triangles  and  hence  each  s  G  <S_oo  is  a  polynomial  of 
degree  <  k  on  R2.  Therefore,  if  T-oo  =  {R2},  then  <S_ ^  a  subspace  of  II*.. 

Suppose  T-oo  =  (A^o)^Tj  and  s  G  S- oo-  Then  each  triangle  A^  can  be  represented  as 
the  union  of  a  sequence  of  nested  triangles  and  hence  s  is  a  polynomial  of  degree  <  k  on  A/,. 
Therefore,  in  this  case,  s  G  <S_ <*,  implies  s  G  G?'(R2)  and  s|Ai  =  Pj\Aj  for  some  Pj  G  11^ , 

j  1  •  V  -  •  •  \  x  • 

Furthermore,  if  s  G  <S_ ^  and  \{x  G  R2  :  |s(x)|  >  t.}\  <  oo  for  some  t  >  0,  then  s  =  const. 
In  particular,  if  s  G  S-oo  H  Lp  (p  <  oo),  then  s  =  0. 
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2.3  Definition  of  B-spaces.  Equivalent  norms.  Interpolation. 

Suppose  T  is  an  LR(or  better)-triangulation  and  <t>  =  $7-  is  a  family  of  differentiable  piece- 
wise  polynomial  basis  functions  over  T  as  described  in  §2.1  2.2.  For  the  characterization 
of  nonlinear  n-term  /^-approximation  from  <4>,  we  need  the  B-spaces  B/ ( $ )  which  we  shall 
introduce  and  explore  in  this  subsection.  In  fact,  the  spaces  //“($)  depend  only  on  the  un¬ 
derlying  ladder  of  spaces  •  •  •  C  <S_i  C  <Sq  C  <S|  C  •  •  •  associated  with  the  bases  but 

as  it  will  be  shown  below  these  spaces  have  atomic  representations  using  <F,  which  justifies 
our  notation. 

We  shall  need  the  B-spaces  -B“(<F)  in  two  cases:  (a)  0  <  p  <  00  and  a  >  0,  or  (b)  p  =  00 
and  a  >  1  (see  Remark  2.14).  In  both  cases,  we  define  r  from  the  identity  1/r  =  a  +  1/p 
( l/oo  :=  0). 

Definition  of  £?“(<!?)  via  local  approximation.  We  define  the  B-space  //“($)  as  the  set 

of  all  functions  /  G  Lr(R2)  such  that 

ll/lk(»)  :=  (  £(|Ar"SA(/)Tr)1/T  <  00,  (2.30) 

a  er 

where  Sa(/)t  is  the  error  of  LT-approximation  of  /  on  OJA  from  Sm  if  A  G  Tm,  (see  (2.22)). 

It  is  readily  seen  that  -B"($)  is  a  linear  space,  ||c/||sa  =  |c|||/||a»  and  \\f  +  g\\ga  < 
II/IIbq  +  II5,IIbq)  with  A  :=  min{r,  1).  Clearly,  see  Theorem  2.8,  if  ||/[|sa  =  0,  then  /  =  0 
a.e.  Therefore,  ||  •  ||#a  is  a  norm  if  r  >  1  and  a  quasi-norm  if  r  <  1. 

We  next  define  other  equivalent  norms  in  _£?“($).  We  define 

JW/)  :=  (  E< \^1,'’SMU)T)1/\  (2.31) 

Acr 


where  we  have  taken  into  account  that  1/r  :=  a  +  1/p.  Thus,  =  ||/||bq($).  More¬ 
over,  we  shall  show  that  ~  ||/||n“(<i>)  if  0  <  7]  <  p  (see  Theorem  2.10). 

Definition  of  norms  in  £?"(<!>)  via  basis  functions  (atomic  decomposition).  For 

/  G  Lr(M2),  we  define 

JV«(/):=  inf  (Ed^riwIW')17'.  <2-32) 


where  the  infimum  is  over  all  representations  of  /  in  the  form  /  =  ^20ee  Cgpg  in  Lr.  (Note 
that  the  existence  of  such  representations  for  each  /  G  Lr  follows  by  Lemma  2.6.)  By 
Theorem  2.9, 


^2i\Eg\  a\\cgpg\\T)r  <  oo  implies 
oee 


oee 


<  00, 

p 


and  hence  Xloee  co^o{^)  converges  absolutely  a.e.  Therefore,  the  specific  type  of  convergence 
that  we  use  in  the  definition  of  N$(f)  above  is  not  essential.  Using  (2.14),  we  have 


AW) 


inf 

/=Xbee  cePs 


(£(lB«l1/f'1/l‘w«IUr)1/' 

oee 
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(2.33) 


~  f_yinff „  (Sllwll ;) ' . 

f-zZoeoCeVe  ^  ' 

Definition  of  norms  in  via  projections.  For  /  G  Ll°c ,  we  set 

Qm,r,(f)  :=  Qrn^if)  ~  Qm-l,r,(f)  e  Sm ,  (2.34) 

where  Qm is  from  (2.21),  and  let  ( be,n{f))oeem  defined  by  the  identity 

QmM)  =  ^  boM)<Po  be,n{f)  ■■=  aoiqmM))’  OeQm.  (2.35) 

oeem 

We  define 

N^Af)  ■■=  ('£,(\Ee\-a\\hlAf)Mrr)',T  (2.36) 

oee 

and,  more  generally  (see  (2.31)), 

W.q, ,(/):=  (5I(|B»r/!'-1/,IIV,(/)^ll„)T)1/T,  0  <T]  <  p.  (2.37) 

6»G0 

By  Lemmas  2.2  2.3,  it  follows  that 

aw/) «  (e  E  tAr/p_1/i'W/)ik<A>r)1/T  (2.38) 

meZ  Ag7T,i 

and,  for  0  <  //,  <  oc. 

w.o.,(/)«  (E(iB»i1/'’_1/'iii'<-.,(/)^ii„)T)1/T«  (Env,(/)^ii;)1/T.  (2.39) 

oee  oee 

We  shall  show  (see  Theorem  2.10  below)  that  all  of  the  above  norms  are  equivalent.  To 
this  end,  we  need  the  following  embedding  theorem. 

Theorem  2.8.  If  f  G  LT(R2)  and  N$tQtr)(f)  <  oo,  0  <  rj  <  p,  then 

f  =  9mM)  =  Y1 W/)^  (2-40) 

mez  oee 

with  the  series  converging  absolutely  a.e.  and  on  Lp,  and 

ll/ll,  <  ell  ElW/KOI  <  c||  ^2\boM)M-)\  <cN^,qM)  (2-41) 

"mez  p  11  oee  p 

with  c  independent  of  f. 

The  proof  of  Theorem  2.8  hinges  on  the  following  more  general  embedding  theorem, 
which  is  a  special  case  of  Theorem  2.5  from  [54], 
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Theorem  2.9.  If  0  <  r  <  p  <  oo,  or  p  =  oo  and  0  <  r  <  1.  then  for  any  sequence  of  real 
numbers  ( co)oeQ  we  have 


oee 


<  c(£llw||;) 7 , 

6>e© 


(2.42) 


with  c  independent  of  (cg)ge@. 

For  completeness,  we  give  the  simple  proof  of  this  theorem  in  the  appendix  (§6). 


Proof  of  Theorem  2.8.  We  introduce  the  following  abbreviated  notation:  Qm  :=  Qm,n(f  )i 
Qm  ■=  qm,n(f ),  bg  :=  bgyT}{f ),  and  N{f)  :=  N^Q^{f).  By  (2.35),  (2.39),  and  Theorem  2.9,  we 


have 


<  c 

p 


£l^(-)l 

eee 


p 


<  cN(f)  <  oo, 


(2.43) 


and  hence  Emez  kmt1)!  <  oo  a.e.  On  the  other  hand,  by  Lemma  2.6,  we  have  ||/  — 
Q,„  |  a:  — >■  0  as  m  — ?•  oo  for  each  A  G  7~.  The  above  two  facts  imply 


OO 

f  —  Qo  =  absolutely  a.e.  on  R2.  (2.44) 

m=l 


We  use  Lemmas  2.1  and  2.2  to  obtain,  for  A  G  77»  ( rn  G  Z), 

IkmlliocfA)  <  ';A  '|</,„|  ma>  <  c|A|“p  £  \\bo<Po\\p  <  <'  A  '..V(/). 

«•  £#DA 


Therefore,  for  a  fixed  A '(=.%,  {v  G  Z), 


E  <  <jv(/)  E  gi 


-i/p 


Aer,  AO  A' 
c,V(/)|A' 


.'I-1*  E  (G'l/GI)1* 

Aer,  AO  A' 

<  c|A,|-1/foY(/)  <  oo, 


where  we  used  (2.9).  We  set 


(2.45) 


Qo  —  E  qm  pointwise  in 


(2.46) 


From  (2.45),  it  follows  that  S'oo  is  well  defined  and  the  series  in  (2.46)  converges  uniformly 
on  every  compact  in  R2.  Evidently,  (2.46)  yields  s0 0  =  Qv  —  oo  dm  for  each  v  G  Z. 

Fix  n  G  Z.  Using  Theorem  2.9,  we  obtain,  for  n  <  n, 


inf  H-Soo  -  .s||p  < 
SE<Sn 


< 


C(  £  I  b>l  r"(J  |  p)  ^0 

fleU™=_oo  ©m 


as  n  — >■  — oo, 
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where  we  used  that  {J^eee  \\bo^o\\Tp)1/r  N(f)  <  oo.  Therefore,  .s-:0C  G  Sn  for  every  n  G  Z 
and  hence  Soo  G  flnez^  =  £-00  • 

Identities  (2.44)  and  (2.46)  yield 

f  ~  soo  =  Yl  =  W/)A?  absolutely  a.e.  (2.47) 

m£Z  0£& 

and  hence,  using  (2.43), 

||/  -  Soollp  <  ell  <  4  IW/)^(-)l  <  cN9,Q,n{f)  <  OO-  (2-48) 

"rnez  p  "see  p 

Since  /  G  Lr  and  /  —  Soo  G  it  readily  follows  that,  for  t  >  0, 

|{z  :  |soo(a:)|  >  £}|  <  |{a;  :  |/(rr)|  >  £/2}|  +  |{s  :  |/(rc)  -  s<x>(£)|  >  t/2}| 

<  (V2)-tii/ii;  +  (t/2)-^n/  -  Soo||j(  <  oo, 

which  implies  .s-:0C  =  0  (see  the  end  of  §2.2).  From  this,  (2.47),  and  (2.48),  we  infer  (2.40) 
and  (2.41).  The  proof  is  complete.  ■ 

Theorem  2.10.  The  norms  ||  •  ||b«($),  (0  <  rj  <  p),  N^(-),  and  N^^q^(-)  (0  <  r/  < 

p),  defined  in  (2.30)  (2.32)  and  (2.37),  are  equivalent  with  constants  of  equivalence  depending 
only  on  p,  a,  rj,  and  the  parameters  of  T  and  <E>. 

Proof.  Theorem  2.8  readily  implies 

AM/)  <  ,'Q^if),  0  <  rj  <  p ,  (2.49) 

if  N^if)  <  °°- 

Suppose  AMs,jj(/)  <  oo.  For  each  A  G  Tm  ( rn  G  Z),  we  have,  by  (2.34)  and  Lemma  2.5, 

|  f/m,/;(./')  |  /.^{  '  }  T  C\\f  ~  Qm,ri  ||l^(A)  T  c\\f  —  Qm—t,ri  ||i^(A)  T  cSa(/)jj  +  cSa«(/)jji 
where  A®  D  A,  A*  G  Tm-it  is  the  only  parent  of  A.  These  estimates  readily  imply 

AW/)  —  Ar$, s,r?(/),  0  <  ij  <  p.  (2.50) 

It  remains  to  prove  that 

Af*,sfl(/)  <  AM//  0  <  r)  <  p,  (2.51) 

provided  AM/)  <  oo.  Evidently,  (2.49)  (2.51)  imply  the  desired  equivalence  of  norms. 

Notice  first  that,  by  Holder’s  inequality,  A'MM/)  <  N$$^(f)  if  0  <  //  <  rj,  and  hence 
it  suffices  to  prove  (2.51)  only  for  r  <  rj  <  p. 

Suppose  /  G  Lr  and  0  <  N$(f)  <  oo.  Then  it  follows  by  the  definition  of  AM/)  that 
there  exists  a  sequence  ( co)oee  such  that 

f  =  ^2cepe  in  Lr  (2.52) 

oee 
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/  \  VT 

and  (  Soe0(l^l_Q|lcoV?6'l|r)T j  <  2A /$(/).  Theorem  2.9  implies  that  in  (2.52)  we  have 
absolute  convergence  a.e.  We  next  estimate 

AW/)  :=  (  £[lAl1/P"1/,,S^a]T)1/T’  (2-53) 

a  er 


using  that  S A{g)v  =  0  if  g  G  Sm  and  A  G  Tm,  and  SA(g)„  <  in  general.  We 

denote  fj  :  =  C<?VW  Fix  A'  G  T  and  assume  that  A7  G  Tm  ( rn  G  Z).  We  have,  using 

Theorem  2.9  (r  <  rj  <  oo)  and  (2.14), 


M/);  = 


< 


< 


*>A' 


oo 

(  E 


< 


j  =m+ 1 


E  fj 

j=m+ 1 


<  c  ^  Hwll) 


oe&,E0cn2t 


c  Y  \E9\ril/T>-1/T)\\c9ip9\\TT. 
oee,Ee  cn% 


Substituting  this  in  (2.53),  we  obtain 

N^UV  <  cY\A'\T(1/p~1/v)  Y  \Ee\T(1/r,-1/T)\\ceMTr 

A'er  ee&,Egcti%t 

=  A  E  (is,i/iA'i^-i/rt(|£#|-«||C,^||r)T 

A'er  ee&,Egcn% 

<  cj:(\E,n\cM\u  E  (|£»i/ia'|)'(i/’>-i/->, 

0e©  a  'eT,n%DEg 


where  we  once  switched  the  order  of  summation.  By  (2.1)  (2.2), 

#{A7  G  Tv  :  Sl%  D  A}  <  c(7V0,  4-  i/eZ,  £  G  0, 

and  l-Efll  <  c/F|A7|  if  A*?  C  with  A7  G  7^  and  6*  G  0m+j  (m  G  Z,  j  >  0).  Using  these, 
we  obtain 

OO 

Y  (|A|/|A7|)r(1/7?-1/p)  <  cY^~1/p)  <  c  <  oo. 
a  'er,n%DEe  j= ° 


Therefore,  N^>,s,r,{f)T  <  cJ2eee (\Ee\  a|lc^e||r)T  <  cN^{f)T  which  yields  (2.51).  ■ 

The  following  embedding  result  is  quite  obvious. 

Theorem  2.11.  For  0  <  «o  <  ot\  and  tj  :=  (ctj  +  1/p)-1,  j  =  0,  l,  we  have  the  continuous 
embedding 

B“/($)c  £«?($),  (2.54) 

i.e.,  if  f  G  ££($),  then  f  G  £“"($)  and  ||/||B«o(4)  <  c||/||B«iW 
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Proof.  By  Theorem  2.8,  if  /  G  then  /  G  LTl  n  Lp  C  LTo.  Fix  0  <  r)  <  p.  Then  by 

(2.39),  we  have 

ll/lls^^)  ~  (XI  IIW/HII?)  /j’  J  =  O’1’ 

3  oee 

and  the  theorem  follows  since  t\  <  tq.  ■ 

Interpolation  of  B-spaces.  We  first  recall  some  basic  definitions  from  the  real  interpola¬ 
tion  method.  We  refer  the  reader  to  [3]  and  [4]  as  general  references  for  interpolation  theory. 
For  a  pair  of  quasi-normed  spaces  X(].  Ah .  embedded  in  a  Hausdorff  space,  the  space  Ah  +  Ah 
is  defined  as  the  collection  of  all  functions  /  that  can  be  represented  as  /o  +  fi  with  /0  G  Ah 
and  fi  G  A, .  The  quasi-norm  in  Ah)  +  Ah  is  defined  by 

ll/IUo+Xt  :=  inf  ||/o||xo  +  H/i|Uv 

J=Jo+h 

Peetre’s  A- functional  is  defined  for  each  /  G  Ah  +  Ah  and  t  >  0  by 

:=  A'(/,£;Ah,Ah)  :=  inf  ||/o||x0  +  ^11/ilUv  (2.55) 

J=lo+h 

The  real  interpolation  space  (X0,  Ah  )\.q  with  0  <  A  <  1  and  0  <  q  <  oo  is  defined  as  the  set 
of  all  /  G  Ah  +  Ah  such  that 

a 00  dt\  1/q 

(t~XK(f,t))q—j  < oo 

with  the  Lq- norm  replaced  by  the  sup-norm  if  q  =  oo. 

It  is  easily  seen  that  if  Ah  C  Ah  (Ah  continuously  embedded  in  Ah),  then  K (/,  t )  ll/llxo 
for  /  G  Ah  and  t  >  1  and,  consequently, 

ll/ll(A-„.A'.),„  «  ll/lk  +  (2.56) 

v=0 

Theorem  2.12.  Suppose  0  <  p  <  oo  and  >  0  or  p  =  oo  and  ao,a\  >  1.  Let 

Tj  :=  (aj  +  1/p)-1,  j  =  0, 1.  Then 

(■B“(*),B"(*))a.t  =  S?(«)  (-«7) 

with  equivalent  norms,  provided  a  =  (1  —  A)o0  +  Aaq  w‘ith  0  <  A  <  1  and  r  :=  (a  +  1/p)-1. 

Proof.  We  shall  use  some  ideas  from  [32],  We  may  assume  that  «o  <  «i-  We  denote  briefly 
Ba  :=  £?“($)  and  Bni  :=  B^d?),  .i  =  0,1.  Furthermore,  we  denote  by  1/  the  space  of  all 
sequences  a  =  (fphee©  of  real  numbers  such  that 

Mk:=(£l“»r)I/,<°°- 

oee 

We  shall  utilize  the  following  well-known  interpolation  result  (see,  e.g.,  [3]): 

(fro,  fTl ) a,t  =  fr,  where  A  =  WA  +  A  with  0  <  A  <  1.  (2.58) 
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We  fix  0  <  rj  <  p.  Then  we  normalize  the  basis  functions  from  <f>  in  Lp,  that  is,  \\^>e\\p  =  1 
(we  use  the  same  notation  for  the  normalized  basis  functions).  We  also  renormalize  the  dual 
functionals  \g  in  the  definition  of  Qm  in  (2.18)  accordingly. 

We  denote  by  b (/)  =  ( M/))fle0  the  sequence  of  numbers  defined  by  (see  (2.34)  (2.35)) 

VmM)  ='■  me  Z  (H^ellp  =  1). 

oe&m 

By  Theorem  2.8,  Theorem  2.10,  and  (2.39),  if  /  G  Bai  ( j  =  0, 1),  then 

/  —  ^ bg(f)(p6  and  \\f\\Baj  ~  ||b(/)||rT. ,  (2.59) 

oee 

and  similarly  for  /  G  Ba. 

The  theorem  will  follow  by  (2.58)  and  the  following  lemma. 

Lemma  2.13.  For  f  G  Ba°  +  Bai  =  Ba°  (a0  <  ot\),  we  have 

K(fJ,]BaFB^)^K(b(f)JrJT0JT1),  t>  0.  (2.60) 

Proof.  We  first  prove  that 

K(f,tBa\B^)<cK(b(f)JrJT0JT1),  t  >  0.  (2.61) 

Indeed,  let  a  =  (a<?)tfe©  £  f’n-  Then  a  G  tTo  (to  >  Ti)  and  since  b(/)  G  tTo  (/  G  Ba° ),  we 

have  b(/)  —  a  G  iT0.  We  define  g  :=  YheeeasFe-  Then  by  Theorem  2.9,  g  is  well  defined, 
and  hence 

./'  -  <1  =  J](M/)  ~  ae)(pe. 

eee 

By  (2.33)  and  Theorem  2.10,  we  infer 

<  clla||rT1  and  \\f  -  g\\B*0  <  c||b(/)  -  a||^T(j. 

Since  a  G  fT1  is  arbitrary,  the  last  two  estimates  give  (2.61). 

We  next  prove  that 

K(b(f)JrJT0JT1)<cK(fJr,Ba\B^),  t  >  0.  (2.62) 

Suppose  g  G  Bai .  Then  by  Theorem  2.11,  g  G  Ba°  (a0  <  cti),  and  hence  /  —  g  G  Ba° .  We 
shall  show  that  there  exists  a  sequence  b (g)  =  ( bo{g))oee  £  f’n  such  that 

g  =  ^be{g)ipe  with  \\g\\B°i  ~  ||b(/)||*Ti  (2.63) 

eee 


and 


f  ~  9  =  ^2(be{f)  -be(g))ipo  with  \\f  -  g\\B*0  «  ||b(/)  -  b(^)||^TQ.  (2.64) 

eee 
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Clearly,  estimate  (2.62)  follows  by  (2.63)  (2.64). 

Notice  that  if  77  >  1,  then  b(-)  can  be  realized  as  a  linear  operator,  and  hence  b(f  —  g)  = 
b (/)  —  b(g).  Therefore,  (2.63)  (2.64)  are  immediate  from  g  E  Bai  and  f  —  g  E  Ba°. 

Suppose  Tj  <  1.  For  A  ET ,  we  let  PA{f)  :=  Pa^U)  £  fl*.  be  the  polynomial  from  the 
definition  of  pm.ri(.f)  in  (2.20)  (PA(f)  is  not  unique).  Thus  PaW  g  n*  is  such  that 

11/  -  -Pa(/)[|l,(a)  <  cEk(f,  A)v.  (2.65) 

We  shall  next  show  that  for  each  A  E  T  there  exists  a  polynomial  P/\{g)  E  IT^  such  that 

II 9  ~  PA(g)\\L,(A)  <  cEk(g ,  A)„  (2.66) 

and 

11/  -  9  -  (PaU)  -  PA(g))hv(A)  <  cEk(f  -  g,  A),.  (2.67) 

We  consider  two  cases. 

Case  1 :  E(f  —  g)  <  E(g ),  where  E(-)  :=  Ek(-,  A)^.  Let  R  E  fl^  be  such  that 

\\f  ~  9  ~  R\\=  E{f  -  g)t  where  ||  •  ||  :=  ||  •  ||^(A).  (2.68) 

We  define  PA{g)  '■=  PaU)  —  R  E  II/..  Then  (2.67)  holds,  by  (2.68).  We  use  (2.65)  and  (2.68) 
to  obtain 

\\g-PA(g)\\  <  c\\f  —  Pa(/)||  +  c\\f  —  g  —  R\\  <  cE(f)  +  cE(f  —  g) 

<  cE(f-g)  +  cE(g)  +  cE(f-g)<cE(g) 


which  gives  (2.66). 

Case  2:  E(g)  <  E(f  —  g).  This  time  we  choose  P&{g)  E  Uk  so  that  \\g~PA{g)\\  =  E(g). 
Similarly  as  above,  one  can  show  that 

Wf-9-(PA(f)-PA(g))\\<cE(f-g). 

Thus  the  existence  of  PA{g)  E  11^.  satisfying  (2.66)  and  (2.67)  is  established. 

Using  the  polynomials  PA{g )  from  above,  we  define,  for  m  E  Z, 

PmM  ■=  1a  '  PUg)  and  pm,n{f  ~  9) '■=  1a  '  (Pa(/)  -  PaU))- 

AeTm  AeTm 

Furthermore,  as  in  (2.21)  and  (2.34),  we  define 

Qm,n(g)  ■ —  QmiPm^ig))  and  Qm,r){9 )  • —  Qm,ri{g)  Qm—  l,r?(fl1)- 

We  define  Qm,v{f  ~  g)  and  qm,v{f  -  g)  in  the  same  way.  Finally,  we  define  b{g)  =  {be{g))oee 
and  b (/  -  g)  =  ( be{f  -  g))0e@  from 

VmM  =:  £  be(gUe  and  qm,r,{f  ~  g)  =■  ^  be{f  ~  g)pe,  m  E  Z. 
oe&m  0eem 

Evidently,  pm,n{f  —  g)  =  pm,ri(f)  ~  Pm,r){g )  and  since  Qm  is  a  linear  operator,  it  follows  that 
b(/  —  g)  =  b (/)  —  b (g).  From  this  and  the  fact  that  PaU)  satisfies  (2.66)  and  (2.67),  using- 
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Theorem  2.8,  Theorem  2.10,  and  (2.39),  we  obtain  that  b(g)  satisfies  (2.63)  and  (2.64),  and 
hence  (2.62)  holds.  This  completes  the  proof  of  the  lemma.  ■ 


By  Lemma  2.13,  (2.58),  and  (2.59)  (with  aj  replaced  by  a ),  we  obtain 


\\f\\(Bao,Bai)XtT  ~  ||b(/)||( 


■^TQ  £t\  )  A  ,T 


l|b(/)lk«  ll/ll 


Ba  ■ 


Thus  the  proof  of  Theorem  2.12  is  complete.  ■ 


Several  remarks  are  in  order. 


Remark  2.14.  (a)  If  p  =  oo,  then  the  B-space  B“($)  (r  :=  1/a)  is  useful  for  our  goals  only 
if  a  >  1.  The  reason  for  this  is  that  Bf(Q)  is  not  embedded  in  C  if  a  <  1.  Indeed,  consider 
the  function  /  :=  Y%  ,  j  1  where  9j  G  @TOj,  m i  <  and  {^  }  are  Courant 

(or  other)  elements  which  overlap  so  that  ||/||oo  ~  Y/  3  1  =  °°-  On  the  other  hand  (see 
(2.33))  |/|B?(*)  <  c{Y%1]~t)1/t  <  oo,  since  r  :=  1/a  >  1. 

(b)  We  introduced  the  B-norms  .Y<i>  )  and  with  0  <  r)  <  p  (see  (2.31)  and 

(2.37))  for  the  following  reason.  As  we  shall  see  in  §3,  normally  a  >  1,  and  hence  r  <  1, 
which  compels  us  to  work  in  Lr  with  r  <  1  that  is  not  a  very  friendly  space.  At  the  same 
time,  if  p  >  1  we  can  choose  1  <  rj  <  p  and  work  in  L v  instead. 

(c)  We  also  want  to  explain  why  we  introduce  the  B-spaces  over  locally  regular  (or 
better)  triangulations  but  not  over  more  general  ones.  The  reason  is  that  if  we  relax  the 
main  conditions  (2.1)  (2.2)  in  the  definition  of  LR-triangulations,  then  we  can  hardly  work 
with  the  B-spaces.  In  particular,  the  equivalence  of  the  norms  (see  Theorem  2.10)  fails  to 
exist,  which  makes  it  impossible  to  prove  all  the  results  from  §3. 


General  B-spaces.  Given  an  LR(or  better)-triangulation  T  and  a  family  of  basis  functions 
T  =  $7-  over  T  as  in  §2.2,  we  define  the  more  general  B-space  ($)  =  B^S),  a  >  0, 
0  <  p,  q  <  oo,  as  the  set  of  all  /  E  Lp(R2)  such  that 


:= 


(E[2”“i 

m£%. 


E 

AeT,  2-™  <|A|<2-m+1 


\  Vrq  <j\ 

Y if)*)  ]  ) 


i\  i/g 


<  oo, 


with  the  f9-norm  replaced  by  the  sup-norm  i f  q  oo .  where  Sa (f)p  is  as  above  (see  (2.22)). 
Evidently,  B“($)  =  R“p($).  In  going  further,  the  norms  in  R“($)  from  (2.31),  (2.32),  and 
(2.37)  can  be  generalized  accordingly.  In  the  present  article,  we  do  not  explore  the  B-spaces 
in  such  generality  because  the  space  scale  R“($)  is  sufficient  for  our  goal  of  characterizing 
the  approximation  rates  of  nonlinear  n-terrn  approximation  from  differentiable  piecewise 
polynomials. 

Fat  B-spaces:  The  link  to  Besov  spaces.  Suppose  T  is  an  arbitrary  SLR-triangulation 
of  R2.  The  fat  B-space  B fk(T)  with  k  >  1  and  a,  r  as  in  the  definition  of  Bf(T)  (§2.3)  is 
defined  (see  [38])  as  the  set  of  all  functions  /  G  Lr(R2)  such  that 

ii/iiiotcn  ~  ( E[iAr"^(/.^wr)I/T  ~  ( E[iAr°^(/.sA)r]T)1/T  <  °o, 

S  a er  a er 
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where  Ek(f ,  I2a)t  is  the  error  of  ^-approximation  to  /  on  Q4.  :=  from  and  LOk{f ,  ^a)t 
is  the  local  Lr-modulus  of  smoothness  of  /  on  Oa-  (Recall  that  Ek(f,  I2a)t  (/i  a)t  by 
Whitney’s  theorem  (2.16),  since  T  is  an  SLR-triangulation.)  Furthermore,  other  equivalent 
norms  in  B“fc(T)  as  well  as  more  general  fat  B-spaces  B“^(T)  can  be  defined  as  in  [38]. 

Suppose  that  <f>  =  $7-  is  a  hierarchical  family  of  basis  functions  over  T  as  described  in 
§2.2.  Assuming  IR  C  Sm  C  Skf  ( T )  for  all  rn  G  Z  (that  is,  k  =  k  in  the  notation  of  §2.2), 
we  have  for  /  G  Lr  and  A  G  Tm, 

Ek(f,tfA)T<c  Ek(f,kl A/)t,  which  implies  ||/||b?(*t)  <  c\\f\\^k{T). 

A'f  7m.  V  <>\ 

Therefore,  the  space  B^T)  is  a  good  candidate  to  replace  B®(Q)  in  nonlinear  spline  approx¬ 
imation,  but  this  is  only  possible  if  0  <  a  <  a$  for  some  «o  <  00,  which  we  do  not  compute 
here.  The  problem  with  the  space  (T)  is  that  ||^||ip*(r)  <  00  only  for  0  <  a  <  a0-  (See 
Theorem  2.15  in  the  case  of  regular  triangulations.)  Therefore,  the  basic  norm  equivalence 
results  (Theorem  2.10)  hold  only  for  a  restricted  range  of  a.  Thus,  B i“fc(T)  is  simply  not 
the  “right”  space  for  the  specific  problem  at  hand  if  a  >  eto-  It  is  too  “fat”.  However,  the 
spaces  B ®k(T)  are  still  noteworthy  since  they  are  less  sensitive  to  small  perturbations  of  the 
triangulation  T  and  are  technically  easier.  We  believe  that  a  situation  will  present  itself 
when  they  will  be  the  “right”  spaces. 

Comparison  between  regular  B-spaces  and  Besov  spaces.  We  begin  by  recalling  the 
definition  of  the  classical  Besov  space  by  moduli  of  smoothness.  So,  the  space  Bsq(Lp)  :  = 
Bsq(Lp(R 2)),  s  >  0,  1  <  p,  q  <  00,  is  defined  as  the  set  of  all  functions  /  G  LP(R 2)  such  that 

a  00  dt  \  1/ ^ 

{t~sujk{f,t)p)qjJ  <00  (2.69) 

{\\f\\p  is  usually  added  to  the  right-hand  side  above),  where  k  :=  [s]  +  1  and  iOk{fB)p  is 
the  A>th  modulus  of  smoothness  of  /  in  Lp(R 2),  i.e. ,  u)k{f,t)p  :=  sup|/; |<#  ||A|(/,  - ) | |p.  It  is 
well-known  that  whenever  1  <  p  <  00,  if  in  (2.69)  k  is  replaced  by  any  other  A:  >  s,  then 
the  resulting  space  would  be  the  same  with  an  equivalent  norm.  However,  the  situation  is 
totally  different  when  p  <  1  and  this  is  a  reason  for  introducing  A:  as  a  parameter  of  the 
Besov  spaces  in  the  following. 

As  elsewhere,  let  us  assume  that  0  <  p  <  00  and  a  >  0  or  p  =  00  and  a  >  1,  and  in 
both  cases  1/t  \=  a  +  1/p.  Let  k  >  1.  We  define  the  space  B^a'k(Lr)  as  the  Besov  space 
B^a(LT)  (see  (2.69)),  where  A:  and  a  are  independent  of  each  other.  These  are  the  spaces 
that  naturally  occur  in  nonlinear  spline  approximation  (see  [53]). 

Suppose  that  T*  is  a  regular  triangulation  of  R2  (see  §2.1).  Then  as  shown  in  [38], 
B “fc(T*)  =  B^a'k(LT)  with  equivalent  norms.  (Notice  that  the  smoothness  parameters  of 
B-spaces  and  Besov  spaces  are  normalized  differently,  a  corresponds  to  2a.) 

Let  us  now  assume  that  $7--  =  {<po}  is  a  family  of  basis  functions  over  T*  as  in  §2.2  such 
that  Uk  C  Sm  C  Skf  ( m  G  Z),  where  r  >  0  and  k  >  r.  As  we  mentioned  above,  the  fat 
B-space  B“fc(T*)  and  hence  the  Besov  space  B^a'k(LT)  is  a  good  candidate  to  replace  the 
B-space  R"($ 7-.)  in  nonlinear  n-term  approximation  from  $7-*.  We  next  spell  out  the  exact 
conditions  for  equivalence  of  the  corresponding  norms. 
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(2.70) 


Theorem  2.15.  Under  the  above  assumptions,  if  0  <  a  <  r  +  1  +  1/p,  then 

B2a'k{LT)  =  £“($r.) 


with  equivalent  norms.  Furthermore,  if  a  single  basis  function  <. pg  G  $7-*  does  not  belong  to 
Cr+1,  then  the  equivalence  is  no  longer  true  when  a  >  r  +  1  +  1/p.  More  precisely,  for  such 
pg  and  a,  \  p„\  a-,  ^  =  00  while  \\pg\\B^(<s>Tp  ~  ||^IIp- 

Proof.  As  we  mentioned  before,  I  ,/'|  />•■•: <I>T.)  <  c||/||H?*(7-.)  for  /  G  B and  also  we 
have  ||/||sp*(r-)  ~  ||/[|B2“-fc(L  )«  exactly  as  in  Theorem  2.25  from  [38].  Therefore, 

II/IIb?(*t*)  ^  c\\f\\B^k(LT)  for  /  e  Bfa’k{LT). 

The  proof  of  the  reverse  estimate  follows  in  the  footsteps  of  the  proof  of  Theorem  2.28 
from  [38]  and  we  shall  indicate  only  the  differences.  Using  the  conditions  on  $7--  and  the  fact 
that  T*  is  regular,  one  can  show  by  straightforward  calculations  that,  for  each  6  G  0(T*), 


f  c|£^(1-(r+1)T)  •  /<  ■■■'■■  if  0  <  t  <  \E0\ll2, 
\  c\Eg\  iff  >  lEgf/2. 


(2.71) 


Moreover,  both  sides  of  (2.71)  are  equivalent  if  pg  does  not  belong  to  Cr+1.  In  going  further, 
one  uses  (2.71)  exactly  as  in  [38]  to  complete  the  proof  of  the  theorem.  ■ 


Remark  2.16.  An  interesting  situation  occurs  when  p  =  00  and  r  =  0.  Then  there  is  no  a 
for  which  (2.70)  holds.  This  is  the  case  when  T7--  is  the  set  of  all  Courant  elements  generated 
by  T*  (a  regular  triangulation). 

Comparison  between  different  B-spaces  and  Besov  spaces.  Suppose  $7-  is  a  family  of 
basis  functions  associated  with  an  SLR-triangulation  T  which  allows  arbitrarily  sharp  angles. 
Then  some  extremely  “skinny”  basis  functions  pg  G  $7-  (with  elongated  level  curves)  will 
occur.  It  is  easily  seen  that  such  functions  have  huge  Besov  norms  (see  [38])  compared  to  their 
Lp-norms  as  well  as  their  R(<f>7-)-norms  (see  Theorem  3.2  below)  for  any  smoothness  a  >  0. 
Therefore,  the  B-spaces  for  such  a  triangulation  are  essentially  different  from  Besov  spaces. 
The  situation  is  quite  similar  when  comparing  two  B-spaces  over  different  triangulations. 
Therefore,  the  B-spaces  change  substantially  with  the  triangulations,  thus  making  the  search 
for  the  “right”  triangulation  mentioned  in  the  introduction  a  meaningful  task.  In  contrast 
to  this,  the  standard  Besov  spaces  can  be  used  only  to  characterize  the  approximation  power 
of  piecewise  polynomials  over  regular  triangulations. 

B-spaces  over  compact  domains.  B-spaces  can  be  introduced  on  an  arbitrary  compact 
polygonal  domain  E  C  R2.  A  substantial  difference  would  be  in  assuming  that  each  tri¬ 
angulation  T  of  E  is  of  the  form  T  =  [J^U0  Tm,  where  T0  is  an  initial  level  (triangulation 
of  E)  and  71,%, . . .  are  consecutive  refinements  of  To-  This  approach  is  important  for  the 
applications  (see  [39]). 

B-spaces  in  dimensions  d  >  2.  Multilevel  triangulations  and  B-spaces  can  be  introduced 
in  much  the  same  way  in  dimensions  cl  >  2.  Of  course,  then  the  triangles  should  be  replaced 
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by  simplices,  thus  making  some  geometric  argumentation  of  this  section  essentially  more 
involved.  In  particular,  the  property  (e)  of  a  multilevel  triangulation  should  be  extended 
to  all  faces  of  the  simplices  in  Tm,  thus  saying  that  there  are  at  most  N0  simplices  in  Tm 
attached  to  a  particular  face.  The  “no  hanging  vertices”  condition  (d)  should  be  replaced  by 
the  condition  that  each  facet  of  a  simplex  in  Tm  is  a  common  facet  of  exactly  two  simplices 
in  Tm.  The  minimal  angle  condition  appearing  in  the  definition  of  regular  triangulations  and 
in  (2.4)  should  be  replaced  by  the  shape  regularity  condition  that  postulates  the  existence  of 
an  upper  bound  on  the  ratio  of  the  diameter  of  a  simplex  and  the  diameter  of  the  inscribed 
sphere.  In  conditions  (2.1)  (2.3)  the  area  should  be  replaced  by  the  d-dimensional  volume. 

B-spaces  in  dimension  d  =  1.  B-spaces  can  be  introduced  in  the  univariate  case,  but 
none  will  give  anything  new  and  hence  they  are  not  needed.  The  key  fact  is  that,  in  the 
univariate  case,  the  Bernstein  inequality  involving  Besov  spaces  holds  with  no  restrictions 
on  the  smoothness  parameter  a  <  oo  (see  [53]). 

In  a  nut  shell,  the  essence  of  the  spaces  we  considered  in  this  section  is  the  following. 
The  Besov  spaces  are  based  on  local  polynomial  approximation  over  regular  multilevel  tri¬ 
angulations,  which  is  explicitly  shown  in  [38].  When  the  regular  triangulations  are  replaced 
by  SLR-triangulations,  then  the  Besov  spaces  become  fat  B-spaces,  which  further  evolve  to 
B-spaces  when  the  local  polynomial  approximation  is  replaced  by  local  spline  approximation. 

The  B-spaces  are  closely  related  to  certain  anisotropic  maximal  functions,  nonclassical 
differentiability,  and  other  problems,  which  are  beyond  the  scope  of  this  article. 

3  Nonlinear  n-term  spline  approximation 

In  this  section,  we  assume  that  T  is  a  locally  regular  (or  better)  triangulation  of  R2.  Also, 
we  assume  that  T  =  $7-  is  a  hierarchical  family  of  basis  functions  over  T  (see  §2.2).  Notice 
that  T  is  not  a  basis;  <E>  is  redundant.  We  consider  nonlinear  n-term  approximation  from 
$  in  Lp(R2)  (0  <  p  <  00),  where  we  identify  L^R2)  as  Co(R2).  We  let  En(<l>)  denote  the 
nonlinear  set  consisting  of  all  splines  -s  of  the  form 

S=£  agipg, 
oeM 

where  M.  C  O(T),  #A4  <  n,  and  M.  may  vary  with  s.  We  denote  by  an(f ,  <f>)p  the  error  of 
Lp-approxi mati on  to  /  €  Lp( R2)  from  En(<f>): 

<7n(/,$)„:=  lnL  11/ “ -sIIp- 

s£E„($) 

Our  goal  is  to  characterize  the  approximation  spaces  generated  by  nonlinear  n-term  ap¬ 
proximation  from  T .  To  this  end  we  next  prove  a  pair  of  companion  Jackson  and  Bernstein 
estimates.  We  shall  utilize  the  B-spaces  £?"($)  introduced  in  §2.  We  assume  that  0  <  p  <  00 
and  n  >  0.  or  p  oc  and  a  >  1.  In  both  cases,  1/r  :=  a  +  1/p  ( l/oo  :=  0). 

Theorem  3.1.  [Jackson  estimate]  If  f  G  Bf(Q),  then 

Vnif,  ®)p  <  cn“a||/||B?(#j  (3.1) 

with  e  independent  of  f  and  n. 
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In  the  case  0  <  p  <  oo,  this  theorem  follows  by  the  general  Theorem  3.4  from  [38],  in 
view  of  the  results  of  §  2.  For  completeness,  we  shall  give  its  short  proof  in  the  appendix. 
The  proof  when  p  =  oo  can  be  carried  out  as  the  proof  of  Theorem  4.1  from  [39]  but  is  a 
little  longer  and  we  shall  skip  it. 

Theorem  3.2.  [Bernstein  estimate]  If  s  E  En(<E>);  then 

INI-B^#)  <  cn“||s||p  (3.2) 

with  c  independent  of  s  and  n. 

The  proof  of  this  (vital  for  our  development)  theorem  utilizes  the  ideas  of  the  proofs  of 
Theorem  3.6  from  [38]  (0  <  p  <  oo)  and  Theorem  4.2  from  [39]  (p  =  oo)  but  is  not  identical 
to  them.  We  shall  give  it  in  the  appendix. 

For  a  fixed  T  and  <E>  :=  $7-,  we  set  K(f\  t)  :=  K (/,  f;  Lp ,  B“($))  ( Lp  :=  Co  if  p  =  00),  see 
(2.55).  The  Jackson  and  Bernstein  estimates  from  Theorems  3.1  and  3.2  imply  in  a  standard 
way  (see,  e.g.,  [55])  the  following  direct  and  inverse  estimates:  For  any  a  >  0,  if  f  E  Lp , 
then 

an{f^)p<cK{f,n-a)  (3.3) 

and 

A-(/,rr“)  <  cn-“  (||/||p  +  [£  bi/V„(/.  $),)"] 1/").  (3.4) 

v=l  U 

where  //,  :=  min{p,  1}  and  c  is  independent  of  /  and  n. 

An  immediate  consequence  of  (3.3)  and  (3.4)  is  that  crn(f,  $)p  =  0(n-7),  0  <  7  <  a,  if 
and  only  if  K(fpn~a )  =  0(n~J).  More  generally,  these  estimates  enable  us  to  characterize 
the  approximation  spaces  generated  by  nonlinear  n-term  approximation  from  T .  We  define 
the  approximation  space  A7  :=  A7($,  Lp),  a  >  0,  0  <  q  <  00,  as  the  set  of  all  functions 
f  E  Lp  such  that 

II/II.43  :=  \\f\\p+  (EtnV-l/.SWd)17’  < 

n=l 

with  the  tq- norm  replaced  by  the  sup-norm  if  q  =  00  as  usual. 

The  direct  and  inverse  estimates  (3.3)  (3.4)  readily  imply  (see,  e.g.,  [55])  the  following 
characterization  of  the  approximation  spaces: 

Theorem  3.3.  If  0  <  7  <  a  and  0  <  q  <  00,  then 

with  equivalent  norms. 

In  one  specific  case  the  interpolation  spaces  can  be  identified  as  B-spaces. 

Theorem  3.4.  Suppose  0  <  p  <  00  and  a  >  0  or  p  =  00  and  a  >  1,  and  let  r  :=  (cx+l/p)-1. 
Then 

A?(c&,  £,)=£“($)  (.3.5) 

with  equivalent  norms. 
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The  following  interpolation  result  is  immediate  from  Theorems  3.3  and  3.4. 

Corollary  3.5.  Suppose  p,  a,  and  r  =:  r(a)  are  as  in  the  hypothesis  of  Theorem  3.  ft,  and 
let  ft  >  a  and  r(ft)  :=  (ft  +  1/p)-1.  Then 

%.*,»  =  (3-6) 

with  equivalent  norms. 


Proof  of  Theorem  3.4.  We  shall  employ  the  idea  of  the  proof  of  Theorem  3.3  in  [30]. 
We  shall  use  abbreviated  notation:  A"  :=  Bft  :=  £?/($),  and  the  alike.  For  any 

ft  >  0,  we  denote  r(ft)  :=  (ft  +  1/p)-1. 

We  first  prove  the  following  continuous  embedding: 

C  with  //.  :=  min{r(/3),  1}.  (3.7) 

Indeed,  suppose  /  G  and  let  sm  G  Em  be  such  that 


Since  crm(f)p  -G  0,  we  have  / 
hence  (//  <  1) 


11/  -  sm\\p  <  2 am(f)p. 

Si  +  Y^L i(s'2''  “  s2*'-1 )  in  (uniformly  if  p 


<  Si 


r(/3) 


E 

V  =  1 


S2"  —  S2"-i 


r(/3) 


(3.8) 


oo),  and 


(3.9) 


We  apply  the  Bernstein  estimate  from  Theorem  3.2  to  s2^  —  s>- 1  G  E2-'+i  and  use  (3.8)  to 
obtain 

1 1 S2^  —  S2*'-1  Hb/3  E  c2V^\\s2»  —  S2"-i||p  E  c2l>^(<7 2’y(f)p  +  cr2"-1(/)p) 

t(B) 

and  similarly  ||si||b/3  E  c(||/||p  +  (Ji(f)p)-  Using  these  in  (3.9)  implies 

t(/3) 


t(/3) 


<c\\f\\;  +  cJ2(^^(fU<4frAS. 

I  ^ 

13=1 


which  is  (3.7). 

Second,  the  Jackson  estimate  from  Theorem  3.1  gives  the  continuous  embedding 

Ect  C  Ai.  (3,10) 

A  third  important  ingredient  in  this  proof  is  the  fact  that  the  approximation  spaces  A“ 
are  invariant  under  interpolation  (see  [31,  52]):  If  «o,a:i  >  0  and  0  <  qi,q2,q  E  oo,  then 

(A“°,  A^)x,q  =  A“,  where  a  =  (1  —  A)cxo  +  Aon  with  0  <  A  <  1.  (3.11) 

Now  we  choose  cto  and  cxj  so  that  0  <  «o  <  a  <  ay  (cko  :=  1  if  p  =  oo).  Also,  we 
select  0  <  A  <  1  so  that  a  =  (1  —  A)«o  +  Ady.  Furthermore,  we  set  r:l  :=  (aj  +  1/p)-1  and 
fij  :=  min{r?,  1},  j  =  0, 1.  By  Theorem  2.12,  we  have 

K°,B“')x,r  =  B?. 
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We  use  this,  (3.7),  (3.10),  and  (3.11)  to  obtain  the  following  continuous  embeddings: 


(H 


ao 

TO  ' 


X,T 


C  (A”AZh,r  = 


A 


a 

T  ? 


which  give  (3.5).  ■ 


Algorithms.  In  [39],  there  are  three  algorithms  developed  for  n-term  Courant  element 
approximation  in  Lp  (0  <  p  <  oo).  These  algorithms  can  be  immediately  adapted  to  non¬ 
linear  n-term  approximation  from  any  family  of  differentiable  spline  basis  functions  $7-  on 
a  compact  polygonal  domain  i?  C  l2.  It  is  an  integral  part  of  our  program  that  using 
the  machinery  of  the  B-spaees,  Jackson  and  Bernstein  estimates,  interpolation,  etc.  devel¬ 
oped  in  this  article,  we  can  prove  that  these  algorithms  achieve  the  rate  of  the  best  n-term 
approximation.  This  aspect  of  our  theory  will  not  be  elaborated  on  here  (see  [39]). 

Approximation  from  the  libraries  {3? 7-} 7-.  An  important  element  of  our  concept  for 
nonlinear  spline  approximation  is  the  introduction  of  another  level  of  nonlinearity  by  allowing 
the  triangulation  T  to  vary.  For  a  given  SRL(or  LR)-triangulation  T,  let  $7-  be  a  family  of 
spline  basis  functions  like  the  ones  considered  in  §2.2.  Now  without  changing  the  nature  of 
the  basis  elements  from  $7-,  we  let  T  vary  and  obtain  a  collection  (library)  of  basis  families 
{$7-} 7-.  We  denote 

°n(f)p  ■  =  inf  <Xn(/,  $r)p, 

where  the  infimum  is  taken  over  all  SLR-triangulations  T  with  fixed  parameters  and  we  also 
assume  that  the  parameters  of  $7-  are  fixed.  The  following  theorem  is  immediate  from  the 
Jackson  estimate  in  Theorem  3.1.  We  shall  assume  again  that  0  <  p  <  00  and  a  >  0,  or 
p  =  00  and  a  >  1,  and  in  both  cases,  1/r  :=  a  +  1/p. 

Theorem  3.6.  If  inf t  \\f\\B*($T)  <  oo,  then 

Vn{f)p  <  <11  n  inf  |  ./'I  W»  (<!>t  i 

with  c  depending  only  on  p,  a,  and  the  parameters  of  T  and  $7-. 

The  ultimate  open  problem  here  is  to  characterize  the  approximation  spaces  generated 
by  {( 7n{f)p }  for  a  given  library  of  basis  functions  {<&t-} 7-. 


Global  smoothness  of  functions:  How  to  measure  it?  Here  we  come  to  one  of  the  fun¬ 
damental  questions  in  approximation  theory  (and  not  only  there)  of  how  the  global  smooth¬ 
ness  of  the  functions  should  be  measured. 

In  the  case  of  nonlinear  n-term  /^-approximation  from  a  single  basis  family  $7-,  a  function 
/  should  be  considered  of  smoothness  a  >  0  if  ||/||n“($r)  <  00.  Then  the  rate  of  n- 
term  Lp-approxi mat i on  of  /  from  $7-  is  0(n~a)  (roughly).  If  we  consider  nonlinear  n-term 
approximation  from  a  given  library  of  basis  families  {<Fr}r  (7~  is  allowed  to  vary),  then  a 
function  /  should  naturally  be  considered  of  smoothness  a  >  0  if  inf 7-  ||/||n«(«i>T)  <  00,  which 
means  that  there  exists  a  triangulation  T  :=  7}  such  that  ||/||s“(#T)  <  00 •  Then  the  rate 
of  n-term  Lp-approxi mat i on  of  /  from  the  library  {Trjr  is  0(n~a).  It  is  crystal  clear  to  us 
that  no  single  (super)  space  can  do  the  job  in  this  case.  It  is  an  open  problem  to  develop  an 
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algorithm  for  finding,  for  a  given  function  /,  an  optimal  (or  near  optimal)  triangulation,  i.e., 
a  triangulation  7/  for  which  /  exhibits  maximal  (near  maximal)  smoothness,  using  the  space 
scale  ).  It  is  also  an  open  problem  whether,  for  a  given  function  /  e  Lp,  there  exists 

a  single  triangulation  7}  such  that,  for  all  n  >  1,  the  n-ternr  L)r approximation  of  /  from 
the  library  {<fy}r  can  be  realized  by  n-ternr  approximation  from  $Tf  and,  consequently, 
characterized  by  the  B-spaces  via  interpolation. 

Another  important  related  issue  for  discussion  is  the  smoothness  of  the  approximating 
tool  $7-  :=  {<p>o}  (T  fixed).  Clearly,  in  nonlinear  approximation,  there  is  no  saturation, 
which  means  that  the  corresponding  approximation  spaces  AJ  are  nontrivial  for  all  0  <  7  < 
00.  Therefore,  the  smoothness  spaces  to  be  used  should  naturally  be  designed  so  that  the 
functions  {ipg}  are  infinitely  smooth  with  respect  to  these  spaces.  This  has  been  one  of 
the  guiding  principles  to  us  in  constructing  the  B-spaces.  Thus  each  basis  function  p>g  G  $ 
is  infinitely  smooth  with  respect  to  the  scale  of  B-spaces  £?“($),  which  is  reflected  in  the 
fact  that  ||^e||Ba($)  <  c||<£0||p  for  0  <  a  <  00  (see  Theorem  3.2).  This  makes  it  possible 
that  in  our  direct,  inverse,  and  characterization  theorems  we  impose  no  restrictions  on  the 
rate  of  approximation  a  <  00  (see  Theorems  3.1  3.4).  Also,  this  explains  the  complete 
success  of  Besov  spaces  in  univariate  nonlinear  spline  approximation  (see  [53])  and  why  Besov 
spaces  are  not  quite  suitable  in  dimensions  cl  >  1.  The  latter  remark  needs  a  few  words  of 
explanation:  First,  by  allowing  triangulations  with  arbitrarily  sharp  angles,  we  allow  very 
“skinny”  basis  functions  with  huge  Besov  norms  compared  to  their  Lp-norms  (see  [38]),  which 
precludes  the  use  of  Besov  spaces  in  such  situations.  Second,  even  when  working  on  regular 
triangulations,  the  use  of  Besov  spaces  is  restricted  by  the  Besov  smoothness  (regularity)  of 
the  basis  functions  (see  Theorem  2.15),  while  B-spaces  impose  no  restrictions  on  the  rates 
of  approximation. 

Spline  wavelets  (prewavelets)  and  frames.  In  the  case  of  uniform  triangulations,  spline 
wavelets  play  an  essential  role  in  practical  algorithms.  It  would  be  desirable  to  have  com¬ 
pactly  supported  wavelet  (prewavelet)  bases  or  frames  generated  by  (differentiable)  spline 
basis  families  $7-  over  LR-  or  SLR-triangulations  T-  To  our  knowledge  there  are  no  con¬ 
structions  of  this  type  available,  as  for  now.  Moreover,  there  is  some  evidence  that  such 
constructions  would  be  too  complicated  and  impractical  for  general  triangulations.  However, 
continuous  spline  prewavelets  on  regular  triangulations  with  uniform  dyadic  refinements  are 
available  from  [21,  34,  58].  (See  also  [47].)  Evidently,  nonlinear  n-term  approximation  from 
compactly  supported  spline  wavelets  or  frames,  generated  by  Courant  elements  or  a  smoother 
spline  basis  family  $7-,  cannot  give  a  better  rate  of  convergence  than  nonlinear  n-term  ap¬ 
proximation  from  $7-.  We  hope  that  efficient  algorithms  for  n-term  approximation  from 
such  families  may  provide  a  substitute  for  wavelet  methods  in  situations  where  the  latter  are 
difficult  to  apply  and,  in  particular,  for  approximation  in  L 00. 

Adaptive  tree  approximation.  This  is  a  method  for  nonlinear  approximation  from  piece- 
wise  polynomials  on  (single  level)  triangular  partitions,  which  has  been  developed  recently 
in  [5,  7].  In  [5],  algorithms  are  developed  which  achieve  the  rate  of  the  best  adaptive  tree 
approximation,  while  in  [7]  the  rates  of  approximation  are  related  to  the  smoothness  of 
the  functions  in  terms  of  Besov  spaces.  There  are  substantial  distinctions  between  this  ap¬ 
proach  and  the  one  in  the  present  article.  Namely,  the  approximation  schemes  from  [5,  7] 
use  “single  level”  piecewise  polynomials  on  triangulations  which  satisfy  the  minimal  angle 
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condition,  while  here  we  use  multilevel  (multiscale)  piecewise  polynomial  bases  over  triangu¬ 
lations  which  allow  arbitrarily  sharp  angles.  Therefore,  the  notion  of  “best  approximation” 
in  [5,  7]  is  quite  different  from  the  one  used  here.  Substantial  progress  has  been  made  in  [6] 
in  applying  the  adaptive  tree  approximation  method  for  numerical  (finite  element)  solution 
of  PDEs. 


4  Construction  of  differentiable  basis  functions 

In  this  section,  we  give,  for  any  SLR-triangulation,  a  construction  of  differentiable  spline 
basis  in  Sff .  r  >  1,  k  >  4 r  +  1,  satisfying  the  conditions  from  §2.2.  In  general,  we  follow 
the  scheme  of  [22],  however,  appropriate  modifications  in  the  construction  and  in  the  proofs 
have  to  be  made  since  we  do  not  assume  that  the  triangulation  is  regular.  In  particular, 
we  replace  the  standard  normal  derivatives  to  the  edges  by  derivatives  in  affine  invariant 
directions,  see  the  definition  of  D ^a)  below.  Since  our  construction  is  also  applicable  to 
non-nested  triangulations  (see  Remark  4.8),  we  formulate  the  results  here  for  a  fixed  level  Tm 
assuming  only  conditions  (a),  (d)  (f),  and  (2.3)  of  §2.1  and  making  sure  that  the  constants 
in  (2.11)  and  (2.12)  depend  only  on  Ayr,  IVo,  and  d2. 

Nodal  functionals.  As  before,  let  Vrn  and  £m  be  the  sets  of  all  vertices  and  all  edges 
of  Tm, ,  respectively.  We  shall  describe  the  basis  functions  for  Sm  =  Sk,r(Tm),  k  >  4 r  +  1, 
with  the  aid  of  the  so-called  nodal  functionals  defined  on  Sk'r(Tm)-  These  are  certain  linear 
functionals  involving  the  values  of  the  splines  and  their  derivatives  at  specific  points  in  M2. 
The  functional  corresponding  to  the  simple  evaluation  of  the  splines  at  (6l2  will  be  denoted 
bv  <fi. 

Of  particular  interest  as  evaluation  points  are  the  vertices  v  G  Vm,  where  we  also  need  the 
derivative  evaluation  functionals  of  type  Sv Df  with  e  being  any  edge  in  £m  emanating  from 
v,  and  Sv Df  ,  where  ei,e2  are  adjacent  edges  emanating  from  v.  Here  D®v  denotes  the 
derivative  of  s  of  order  a  in  the  direction  of  the  interval  [i>,7],  weighted  with  the  length  of 
[y,  7],  namely, 


D[V,v]S  ((w  VX)DX  +  ( Vy  Vy)Dy)  s, 

V  =  {vx,vy),  r  ( 'by)- 

Note  that,  due  to  this  weighting,  the  corresponding  Markov  inequality  reads  as  follows: 

—  C||/J||ioo[^,«]’  P  ^  (4-1) 

where  c  depends  only  on  k  and  a. 

Let  Ai,A2  G  Tm.  share  an  edge  e.  Since  every  s  G  Sk,r(Tm )  is  continuous,  the  two 
polynomial  patches  and  -s|a2  coincide  along  e.  Therefore,  5vDfs  may  be  computed  for 
any  a  =  0,1,...  as  either  5vDf(s\ aJ  or  $vDf{s  |a2)  with  the  same  result.  Similarly,  let 
ei,e2  G  £m  be  two  edges  of  a  triangle  A  G  Tm  with  a  common  vertex  v.  Then  SyD^D^s 
denotes  the  mixed  derivative  of  s  at  v  in  the  directions  of  e1  and  e2  away  from  v.  If  a  +  fi  <  r, 
this  derivative  is  uniquely  defined.  If  a  +  0  >  r ,  the  result  may  depend  on  the  choice  of  the 
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polynomial  patch  of  s  attached  to  v.  We  follow  the  convention  to  always  take  dvDf  Dfs  :  = 
5VD£  D@  (s |  a),  where  A  is  the  above  triangle  formed  by  ei,e2. 

We  shall  also  need  functionals  evaluating  at  some  points  on  any  edge  e  the  derivatives 
of  the  spline  in  an  affine  invariant  direction  not  parallel  to  e.  Let  e  =  [t>i,t>2]  €  Sm,  and 
let  Ae  =  [vi,V2,  v3\  G  Tm  be  a  triangle  attached  to  e.  Denote  by  p(e,  A)  the  median  of  A 
connecting  the  middle  point  (tq  +  u2)/2  of  e  with  the  third  vertex  y3  of  A.  For  any  point 
(6e,  dcI)ll!e  j(K]  will  denote  the  derivative  at  f  in  the  direction  pointing  into  the  half-plane 
containing  A  parallel  to  p(e,  A),  weighted  with  the  length  of  p(e,  A).  For  each  edge  e  G  £m . 
we  choose  one  of  the  two  triangles  attached  to  e  and  denote  it  by  Ae.  (Note  that  this 
selection  of  Ae  is  not  unique  but  as  will  be  seen  it  will  cause  no  problems  for  the  basis 
construction.) 

Remark  4.1.  For  later  references,  we  note  here  that  any  nodal  functional  rj  :  Sk'r(Tm)  — >■  R 
of  the  above  type  can  be  extended  to  a  linear  functional  rj  :  Sk~1(Tm)  — >■  R  such  that 
fj(s)  =  p(s)  as  long  as  s  G  Sk,r(Tm).  Indeed,  if  the  definition  of  p  involves  5%  for  some 
point  £  G  Ueesme,  then  we  choose  one  of  the  triangles  A  G  Tm.  containing  f  and  use  the 
corresponding  value  of  s|a  or  its  derivatives  at  £  to  define  rj(s )  for  any  s-  G  Sk~l(Tm).  The 
only  restriction  on  the  choice  of  A  is  that  it  must  be  consistent  with  the  above  rules  for 
<w,  and  S^D^a).  Clearly,  the  extension  of  this  type  is  not  unique.  Moreover, 

convex  combinations  of  evaluations  of  the  restrictions  of  s  to  different  triangles  can  also  be 
used. 


Characterization  of  differentiability.  Let  L  be  a  straight  line  dividing  R2  into  two 
half-planes  H.H.  Given  p.  p  G  IIfe ,  we  define  a  piecewise  polynomial  function  s  by  setting- 
sir/  =  A  s\h  =  P-  To  check  whether  s  is  differentiable  across  L,  we  choose  two  points  u,  v  on 
L.  as  well  as  two  points  w,  w  in  the  interiors  of  H  and  It.  respectively.  We  set  A  :=  [u.  v.  w ], 
A  :=  [u,v,w\,  e  :=  [u,v],  p  :=  [u,w\,  p  :=  [u,w\,  6  :=  Aep,  6  :=  Zpe.  The  proof  of  the 
following  lemma  can  be  found  in  [17,  25]. 

Lemma  4.2.  Let  0  <  r  <  k.  Then  s  G  C"'(R2)  if  mid  only  if 

S.Df.Di-p  =  £<-l)»(*)  sin“-»(9  +  d)(^y(^yaSuDVDl-l>p  (4.2) 

P=0 

for  all  a  =  0, . . . ,  r  and  q  =  a, . . . ,  A:  —  1. 

It  is  readily  seen  that  (4.2)  can  be  reformulated  as  follows: 

ot  a—j3 

(*d$di~*p = Et-n*©  (*|A'i)  (4.3) 

4=o 

where  a  :=  sgn  sin(#  +  9)  and  A*  :=  [u,w,w\.  (This  identity  simplifies  in  an  obvious  way 
when  | A* |  =0). 

See  [22]  for  a  discussion  of  the  relationship  between  these  nodal  conditions  of  differentia¬ 
bility  and  the  well-known  Bernstein-Bezier  conditions. 
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Construction  of  basis  splines.  Consider  the  following  set  Mm  of  nodal  functionals  on 

Aft,  ■.=  (  U  Wk)  u  (  U  Aft)  u  (  U  Aft),  I4-4) 

VdzVm  6G^m  AG7m 

where  for  each  A  =  [tii,  t>2,  ^3]  G  Tm, 


-A  •  — 


j*m  +  *2  A 


AC  :=  {>)f  !=  if  :  «  £ 

■  i 


-A 


h 


k  -  1 


k  +  '■>  +  /3  =  k-1,  h,  /■>.  >  r 


}  C  A, 


for  each  edge  e  =  [  ci,  a]  G  £m, 


Je,g 


■C  :=  {nU  ■=  kDl(eAe)  :  9  =  0,  •  •  • ,  r,  f  G  ~e,  J, 

:=  {  +  ?'2^2  ;  7  +  i2  =  k  -  q  -  1,  *!,  ?;2  >  2r  -  (?)  C  e, 

L  K  —  q  —  1  J 


and  for  each  vertex  v  G  Vm, 


2r 


J\fv  :=  I  I  .AAA 

J  v  m  J  v  m  5 

<7=0 


with  A/j^’9,  9  =  0,...,  2r,  being  defined  as  follows.  Let  A^  =  [u,  t a,;+i],  i  =  1 .... .  7VW,  be 
the  triangles  in  7m  attached  to  v  in  counterclockwise  order,  vNv+£  =  /y.  and  let  e*  =  [n,:§|]. 
We  set 


K°  ■■=  W'°  ■■= 


v,0 


Mi 


K«  ■■=  := 


—  Q;  7~)Q; 
^ei+i 


i  =  1, 


W,  a  =  0, ...  ,q  -  1},  q  >  1. 


Note  that  A/^  or  A/"7  might  be  empty  for  some  combinations  of  r,  k,  e.g.,  A/^  =  Mren  =  0 
if  r  =  0,  k  =  2,  or  A/^  =  0ifr  =  l,A;  =  6.  This,  liowever,  does  not  cause  any  problem  for 
the  construction  below. 

In  view  of  (4.2),  the  functionals  in  J\f^q  are  not  linearly  independent  on  Sk'r(Tm)  if  q  >  1. 
Namely,  the  following  conditions  hold  for  all  s  G  Sk'r(Tm),  v  G  VTO,  9  =  1,...,  2 r: 


■«(»)  -  E(-D"©  +  ft)  (iT^)"(fcdRT7r1)  Ac, =  o, 

3=0 

a  =  1, . . . ,  min{r,  9},  i  =  1, . . . ,  Nv, 


(4.5) 


where  0,  :=  ketel+1,  r%$  :=  %;%. 


The  following  key  lemma  is  instrumental  in  constructing  the  basis  functions. 
Lemma  4.3.  There  is  a  unique  spline  s  G  Sk,r(Tm )  such  that 


7a(s) : 

II 

a 

•n>  O 

< 

01 

UJ 

Uy> 

A  G  Tm, 

Vq^s) 

II 

a 

►O  Ct> 
(go 

<iT 

[i] 

UJ 

q  =  0,  •  •  • ,  r,  e  G  £m, 

rjv'°(s) 

=  a®’0, 

v  G  Vm, 

&s) 

II 

s 

7s-  e 

P  ^ 

• ,  K,  0  9  1.  9  =  1,. 

•  •  1  2r,  i)  G  Vm, 
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for  any  given  af,ae  ^av'°  G  R  and  any  afq  G  R  satisfying 

-i  n;) 

5=0 


av'q  - 

ut,a 


±(-w® 


-  n 

-1.9-5  —  u’ 


a  =  1, . . . ,  min{r,  </},  /  =  1, . . . ,  iV„. 

Moreover,  for  each  A  £  Tm, 

IMaIIwa)  <  cSf2r  max  |r?(s)|, 

■  ♦  ijeMn(A) 

where  c  is  a  constant  depending  only  on  k,  and 

JV™(A):=(  U  •AC)U(  U^LIJV, 


A 

m  ' 


«evmnA 


e££)7 

eCA 


(4.7) 


(4.8) 


Proof.  We  first  determine  s|e  for  each  e  =  [17,  112]  G  £m  using  the  fact  that  s|e,  as  a  univariate 
function  on  the  interval  e,  is  a  polynomial  se,o  of  degree  at  most  k  —  1.  Therefore,  570  is 
uniquely  determined  by  the  following  k  Hermite  interpolation  conditions: 


SVl9e,o  =  aVl'°,  0  =  fl1'2’0, 

SVlDfsefi  =  a^o’7,  SV2UJse, 0  =  ftj20’7,  7=1,...,  2 r, 

<Vse,0  =  ftog)  £  e  1— ‘e,0i 


(4.9) 


where  we  assume  that  e  is  the  z-tli  edge  emanating  from  77  and  the  jth  edge  emanating  from 
1-2  ■ 

We  next  determine  :=  {D9^  _^s)|e,  =  Let  Ae  =  [77,%, 'Vs]-  Then 

Dl(eAe)  =  Therefore,  for  7  =  0, ... ,  2r  -  9, 


A  D7s 


£(-d 


7-f 


KDe*‘Dl~^S  = 


£(-d 


Z  n,q—l  1 


•s). 


r=o 


f=0 


Similarly,  since  =  (£>[t,2lt,3]  -  we  have  for  7  =  0,...,  2r  -  q, 

=  £(-i)'2 -'ABy'-oyy 

£=0 

=  £(-i)'2-SJif7A')  +  (-d«2-«77+54 

f=0 

In  addition,  we  have  for  each  £  G  Hew, 


^Se,9  =  hDl(e,0^S  =  rilz(s). 
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Thus,  for  each  q  =  1 ,r,  the  univariate  polynomial  se,q  of  degree  k  —  1  —  q  is  uniquely 
determined  by  the  k  —  q  Hermite  interpolation  conditions 


=  £( 7  =  0,...,  2 r  -  q. 


1=  0 
q-1 


S„DJse,,  =  £(-l)'2 -'o?t+T  +  7  =  0 . 2r  -  q, 

£=0 

^£se,q  =  ^ 

Let  A  =  [t)i,ti2,«3]  G  Tm  be  the  second  triangle  attached  to  e.  We  set 


(4.10) 


a0,£’ 


£  £  HfcOl 

C—(\  \  / 

r, 


e=o 


Je,?i 


where  A*  :=  [fj,  t>3,  f53],  U  =  (tq  +  t>2)/2,  and  a  :=  sgn  sin(i>3?m2  +  73UU2),  where  uvui  denotes 
the  angle  determined  by  u,v,w  with  vertex  at  v.  (It  may  happen  that  |A*|  =  0.)  Since 
A*  C  conv  (Ae  U  A),  we  have 


A*\q~l\A\l\AP\~q  <  5iq+£\A\q\AP\-q  <  Sk2q+i. 


2i  u2 


(4.11) 


We  now  construct  each  polynomial  patch  .s|A,  A  G  Tm ,  of  the  spline  s  as  the  unique  solution 
of  the  following  interpolation  problem: 


4(SU)  =  "c  . 

< 

[1] 

UJ 

hDl(e,  A)(S'U)  = 

£  £ 

q  =  0, . . . ,  r,  e  C  A 

if  Ae  =  A, 

Aiy!,^A)(s  '  1  =  "'if 

£  G  ^e,q, 

q  =  0, . . . ,  r,  e  C  A 

if  Ae^  A, 

A(.s  \)  =  av'°, 

v  G  A, 

aD*.+i{s  )  = 

a  =  0, . . 

. ,  q  I , .  .-.  .q  2?’,  u 

G  A, 

(4.12) 


(i  is  such  that  ej,ei+i  C  A). 


Since  (4.12)  is  a  standard  finite  element  interpolation  scheme  for  bivariate  polynomials  of 
degree  A:  —  1  (see,  e.g.,  [57]  or  Lemma  3.7  in  [25]),  the  polynomial  s|A  is  uniquely  determined. 

We  now  show  that  the  piecewise  polynomial  s  constructed  in  this  way  lies  in  the  space 
Sk'r(Tm ),  be., it  is  r  times  differentiable,  To  this  end  we  consider  any  edge  e  =  [t>i,t>2]  G  £m. 
As  before,  let  Ae  =  [rq,  «2,  a]  and  let  A  =  [f2,  Wr^s]  be  the  second  triangle  attached  to  e, 
and  we  again  assume  that  e  is  the  7tli  edge  #4,*  emanating  from  v\  and  at  the  same  time  the 
jth  edge  (:‘2,j  emanating  from  V2-  Then  we  have 

ei,i-i  =  [ui,t?3],  ehl  =  [ui,u2],  eiii+i  =  [ui,u3], 
e2,j- 1  =  [  A ,  Ai] ,  e2,j  =  [v2 ,  Ui] ,  e2,j+l  =  [  A ,  A] . 
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Obviously,  for  each  q  =  0, . . . ,  r,  Dq^e  Af)(s'|Ae)|e  =  se.q  satisfies  the  interpolation  conditions 
(4.9)  if  q  =  0  or  (4.10)  if  q  >  0.  We  set 


Je,q 


■=  Dq  -  | 

Me,  A )' 


8 


The  desired  differentiability  of  s  will  follow  if  we  show  that 


8e,q  ~  $e,q  ^(e,A)  a)  lej  ll  0.  .  .  .  .  A 

By  (4.12)  we  have 


C  (Sl  Ae  )  — 

C  ( s  A ) 

=  a”1’0, 

KDt:*K,i+M±')  =  <Cq’ 

a  =  0, 

•••,9-1, 

Q  = 

1,.. 

. ,  2 r. 

KDq~ZD«Js\ A)  =  la , 

a  =  0, 

•••,9-1, 

Q  = 

1  • 

. ,  2r. 

which  in  view  of  (4.7)  imply 


(4.13) 


=  snDt;;D-m(si  A. 

a  =  0, . . . ,  min{r,  9},  q  =  0, . . . ,  2 r, 


and  hence 

KD2{§e,q  -  se,q)  =  0,  7  =  0,...,  2r  —  9,  q  =  0,...fr. 

Similarly,  we  get 

&V2  DZ{se,q  -  se,q )  =  0,  7  =  0,...,  2 r  -q,  q  =  0, . . . ,  r. 

In  addition,  a  simple  calculation  relying  on  (4.3)  shows  that 

S^setq  =  aeq  g,  £  G  S e,qi  q  =  0, . . . ,  r, 

so  that  by  (4.12), 

—  ^e,ij)  =  0,  ^  G  q  =  0, . . . ,  r. 

Since  $(,.q  —  se^q  satisfies  homogeneous  interpolation  conditions  of  a  well-posed  Hermite 
scheme,  (4.13)  follows. 

The  uniqueness  of  s  is  clear  from  the  above  proof,  since  s  =  0  if  the  numbers  in  the 
right-hand  side  of  (4.6)  are  all  zeros. 

It  remains  to  prove  (4.8).  Since  se,q  satisfies  the  interpolation  conditions  (4.9)  if  q  =  0  or 
(4.10)  if  9  >  0, 


I We,q 1 1 Loo te)  <  crnax{?j(.s)  :  r]  G  AC  U  AC  U  A/£},  9  =  0,...,  r, 

where  c  depends  only  on  k.  In  view  of  (4.11)  and  Markov  inequality  (4.1),  we  have 

Kgl  <  cS22q\\Se,q\\Loo(e),  Q  =  0,  .  •  •  , 

and  (4.8)  follows  by  the  properties  of  the  interpolation  problem  (4.12).  see  Lemma  3.9  in 
[25].  ■ 
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For  each  v  G  Vm  and  q  =  1, . . . ,  2 r,  we  denote  by  R^q  the  (min{r,  q}Nv  x  gA^)-matrix  of 
differentiability  conditions  (4.5).  Let  the  vectors 

av'q'\  j  =  1, . . . ,  pv,q  :=  qNv  -  rank {R^q), 

form  an  orthonormal  basis  for  the  null  space  of  Rff1 . 

null(^)  :=  {a  G  RqN “  :  Rvmqa  =  0}. 

For  convenience,  we  shall  use  the  double  indices  introduced  in  the  definition  of  J\fffq  also  for 
the  components  of  av,q'q: 


avfaJ,  *  =  1 . Nv,  a  =  0,  1. 

(4,14) 

We  set 

Nv  q-1 

,nv,qj  _  V  V  av'q,jnv,q  7-1  o 

V  ■  /  ,  /  J  ai,a  ’ll, a'  j  ■  ■  ■  1  Pv,qy 

i=  1  a= 0 

(4,15) 

Kn9  :=  { VV '9'3  ■  j  =  !,•••,  Pv,q},  q  =  1,  ■  ■  •  ,  2 r, 

2  r 

V’:=K:°uUW.  t.€V», 

<?  I 

V„  :=  (  (J  v;„)  U  (  U  V‘)  u  (  U  Ae), 

tlGVm  6G^m  AGTm 

and  define  the  set 

Tj,,  —  ^  W"m  } 

of  the  basis  functions  for  Sk'r(Tm )  by  the  duality  condition, 

J  1  if  n  =  rj , 

\  o  if peAfmMv}. 

(4.16) 

To  see  that  the  above  definition  is  correct  we  have  to  check  that  for  each  rj  G  Mm  there 
exists  a  unique  (pv  satisfying  (4.16).  This  follows  by  Lemma  4.3.  Indeed,  since  the  vectors 
av'q'q  are  orthonormal,  we  have 


v.q 

Via  = 


Pv,q 

E 

3=1 


nv'q'Jnv 
ui,a  '/  ’ 


i  =  l,...,Nv,  a  =  0, . . .  ,q  -  1. 


Therefore,  for  a  fixed  //.  the  numbers 


av'q  — 

ui,a 


■= 


i  =  <\  ()..., -W  1. 


satisfy  (4.7),  which  ensures  the  applicability  of  Lemma  4.3. 
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Properties  of  basis  splines.  It  follows  by  Lemma  4.3  that  every  spline  s  G  Sk'r(Tm )  is 
uniquely  determined  by  the  sequence  (rj(s))  e A'-m,  i.e.,s-  has  a  unique  representation 


S  —  'y  '  ClrjLPrj^ 
rteMrn 


av  =  Tj{s)  G 


Furthermore,  (4.8)  immediately  implies 

star('c)  if  rj  G  J\f^  for  a  vertex  v  G  Vm, 
supp  ( pv  C  ^  star(e)  if  rj  G  J\ffn  for  an  edge  e  G  £m , 

if  ij  G  for  a  triangle  A  G  Tm,, 


A 


r— 2r 


||^||loo(k2)  <  c^2  . 

By  using  Markov  inequality  it  is  easy  to  show  that 

1 3 1 1  Loo  (star  (v))  if  V  G  A C  for  a  vertex  v  G  VTO, 

I S 1 1  Loo  (star  (e))  if  V  e  A/£  for  an  edge  e  G  £m, 

l-sIUoo(A)  if  T)  e  AC  for  a  triangle  A  G  Tmi 


M«)l  <  c  < 


(4.17) 


(4.18) 


(4.19) 


with  c  a  constant  depending  only  on  k.  r  and  JV0. 

Thus,  we  showed  that  the  basis  <Lm  =  {c :  r/  G  A fm}  satisfies  all  requirements  of  §2.2 
with  Sm  =  Sk'r(Tm)  and  k  =  k.  (Obviously,  Ifo  C  Sm  and  Sk'r(Tm)  C  Sk'r{Tm+ 1)  if  Tm+i  is 
a  refinement  of  Tm.)  More  precisely,  we  have  the  following  result. 


Theorem  4.4.  Let  r  >  0,  k  >  4r  +  1.  Suppose  that  Tm  satisfies  (a),  (d)  (f),  and  (2.3)  of 
§2.1.  Then  the  basis  functions  ipv  G  Sk'r(Tm)  (rj  G  J\fm)  constructed  above  have  the  following 
properties: 


(a)  For  any  s  G  Sk,r(Tm. )  there  exists  a  unique  sequence  of  real  coefficients  (g^CaA  such 
that 

^  ]  ar/  Tr)  i 


S  = 


rieAfn. 


with  av  =  Tj{s),  rj  G  A fm. 

( b )  For  each  rj  G  A fm  there  is  a  vertex  v  =  vv  G  Vm  such  that 


supp  ipn  C  star(ri)  =:  Ev, 

IIvaIIwsa  =  IKIIw£a  < 

'/(*)  <  A^II-sIIloo^),  S  G  Sk'r{Tm), 
where  M\.  M_,  are  positive  constants  depending  only  on  k,  r,  S2,  and  N0. 

In  particular,  by  the  proof  of  Lemma  2.3,  we  have  the  following  stability  property  of  $TO. 
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Theorem  4.5.  The  basis  $TO  is  Lp- stable  for  all  0  <  p  <  oo,  i.e.,for  any  sequence 

I  ip(K,2)  ~  (  ^2  WavTv\\PLpm)  5 

riEjifm  VOTfm 

where  the  constants  of  equivalence  depend  only  on  p,  A:,  r,  c>2;  and  N0.  In  the  case  p  =  oo  A/ie 
tp-norm  in  the  right-hand  side  is  replaced  by  the  sup -norm. 

The  linear  functionals  Xv  :  Sk~1(Tm)  H  LQO(ET1 )  — >■  M,  rj  G  Am,  with  properties 

\(s|e,  )  =  Ms),  s  e  Sk'r{Tm), 

l\(/)|  <  Afe||/IU„(j!,).  /  e  Sk-'(T,„)\E,  n  LX(E„), 

needed  in  the  definition  of  the  projector  Qm  (see  (2.18))  can  now  be  defined  in  a  constructive 
manner.  Indeed,  we  first  extend  each  functional  rj  G  J\fm  to  a  functional  fj  defined  on 
Sk~l{Tm ),  according  to  Remark  4.1,  and  then  set 

\  '=  h  if  V  e  (  1J  ^rn)  U  (  U  •A/^) 

AgTJji 


Nv  ?-l 


A  •-  av'q,jfiv'q 

Ar>  ■—  na  h,a 


Nv  <?-l 


2=1  a=0 


/  j  /  j  2, a  ii, a 

2=1  a=0  vEV' 


u  < 


By  (2.22),  Qm  can  be  extended  to  the  operator  Qm.p  :  Ll°c  — >■  Sk,r{Tm )  whose  local 
approximation  power  is  described  in  the  following  theorem  (see  Lemma  2.5). 

Theorem  4.6.  Suppose  f  G  Lloc,  0  <  p  <  oo  (/  G  C  if  p  =  oo).  Then 


\\f  -  Qm„p{f)\\Lp(A)  <  eSA(f)p  <  cEk{f,QA)p,  A  G  Tmi 


where  is  the  union  of  all  triangles  in  Tm  that  have  a  common  vertex  with  A,  and 

the  constant  c  depends  only  on  p ,  A:,  r,  52,  and  N0. 


To  show  that  the  assumption  that  the  triangulations  Tm  satisfy  (2.3)  cannot  be  omitted, 
we  consider  the  following  example. 

Example  4.7.  Suppose  Tm  has  an  edge  e  =  [v,  u]  with  two  triangles  A  =  Ae  =  [«,  u,  w]  and 
A  =  [v,  %  w]  attached  to  e,  such  that  u  =  v  +  (2  At.  0),  w  =  v T  (— a,  a),  w  =  v  +  (— a,  —a), 
where  the  positive  numbers  M,  a  depend  on  m.  Evidently,  |conv  (AeU  A)|/|  Ae|  =  2(2M  +  1), 
and  (2.3)  will  be  violated  if  M  grows  unboundedly  with  rn.  while  the  maximal  angle  of  the 
two  triangles  is  37t/4,  thus  allowing  the  maximal  angle  condition  (2.5)  to  hold.  Note  that 
such  configurations  of  triangles  are  possible  for  a  sequence  of  levels  of  an  LR-triangulation  T 
with  the  corresponding  M’s  tending  to  infinity,  see  §2.1  of  [38].  Choosing  A:  =  6  and  r  =  1, 
we  consider  the  basis  functions  tpri  G  S6'1(Tm, ),  rj  G  J\fm,  constructed  according  to  the  above 
algorithm.  We  next  show  that  the  basis  =  {pv  :  rj  G  Am.}  is  instable ,  i.e. ,  Theorem  4.5 
does  not  hold  for  it.  (Therefore,  neither  $TO,  nor  a  renorming  of  it  satisfies  the  requirements 
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of  Section  2.2.)  More  precisely,  we  show  that  the  constant  function  lK2(a;)  e1,igR2,  does 
not  have  an  Loo-stable  expansion  with  respect  to  <Lm.  We  have 

ll^-K2  IUoo(K2)  =  1,  liR2  = 

Now  choose  rj  =  rf'°  =  8%.  G  Since  77(1^2)  =  1,  the  instability  of  (l>rn  will  follow  if  we 

show  that  ||Lj?IUoo(k2)  is  unbounded  as  M  — >■  00.  By  (4.12), 

y  1  =  au 

=  —  2  —  S^Dlsefi  ~ 

where  £  =  (v  +  u)/ 2,  A*  =  [£,«;,«;],  se,o  =  A?U.  se,i  =  (^' ;, .  ■  jW,)  Obviously,  |A|/|Ae|  = 

1,  and 

I  A*  I/I  Ae|  =  ( |  conv  (Ae  U  A)|  -  /|Ae|  =  2M+1  +  1. 

The  univariate  polynomial  se,o  of  degree  5  is  determined  by  the  Hermite  interpolation  con¬ 
ditions  (4.9)  that  take  in  our  case  the  form 

fivSe,  0  1?  A'^r.o  A  I)v  Sf.j)  8uDesefi  SvDesefi  8VD  ese^  0. 

An  elementary  computation  shows  that  8^D\se$  =  —15/8.  By  (4.10),  we  immediately  get 
<^se,i  =  aeu  =  =  0.  Thus, 

5tDl(e, A)KIa)  =  -j(2"  '  '  +  1)  ^  00  as  M  >  00. 

In  view  of  Markov  inequality,  |  aJ  /  X.(T)  —  c^-D1.  -  (<^jjjT)|,  and  we  get  the  desired  un- 
boundedness  of  HA7IU00A2)  f°r  sufficiently  large  M. 

Remark  4.8.  It  is  clear  that  Theorems  4.4  4.6  are  valid  for  any  sequence  of  levels  Tm 
satisfying  the  hypotheses  of  Theorem  4.4,  i.e.,  nestedness  and  other  additional  assumptions 
on  {Tm}  stated  in  §  2.1  are  not  needed  for  these  results. 

Remark  4.9.  It  is  an  important  property  of  the  basis  functions  tpn  constructed  above  that 
they  are  invariant  under  affine  transforms.  More  precisely,  let  Tm  satisfy  the  hypotheses 
of  Theorem  4.4,  and  let  A  :  M2  — >■  M2  be  an  affine  transform.  We  set  A(Tm)  =  {A(A)  : 
A  G  Tm},  and  Ar)(s)  :  =  8A(v)D®(e^DaAf.e^s  for  each  nodal  functional  7]  of  the  form  7] (s)  = 
8vD“D%  s,  and  extend  the  operator  A  linearly  to  the  linear  combinations  of  the  nodal 
functionals  such  as  those  occurring  in  (4.15).  Then,  clearly,  the  sets  of  nodal  functionals 
Mm  and  N,'}  defined  by  (4.4)  for  Tm  and  A(Tm ),  respectively,  satisfy  N,'}  =  {At)  :  rj  G  J\fm}- 
(We  used  here,  in  particular,  the  fact  that  n(A(e),  A( Ae))  =  A{fi(e,  Ae)).)  Moreover,  since 
the  matrices  of  the  differentiability  conditions  (4.5)  are  affine  invariant  (see  (4.3)), 
we  also  have  J\f^  =  {Arj  :  r)  G  Mm }  for  the  appropriate  sets  Afm,Afm  defined  as  in  the 
construction  above,  provided  we  choose  the  same  orthonormal  vectors  (4.14)  in  both  cases. 
Let  now  =  {a,  :  V  e  J\fm }  C~Sk'r{Tm)  and  =  {ipAri  :  7]  G  J\fm}  C  Sk’r(A{Tm ))  be 
the  spline  bases  dual  to  J\fm  and  A/m' ,  respectively.  Since  r)  G  Mm .  obviously  satisfy 

the  same  duality  relations,  we  conclude  that  A'?/  =  A?  (4b),  r)  G  J\fm,  which  is  the  desired 
affine  invariance. 
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Remark  4.10.  Our  construction  is  extendable  to  the  spaces  Sk'r(Tm, ),  k  >  r 2d  +  1,  in 
dimensions  cl  >  2.  To  this  end  the  algorithm  given  in  [22]  should  be  extended  to  SLR- 
triangulations  in  Rd.  In  particular,  the  orthogonal  directions  of  derivatives  used  in  [22] 
should  be  replaced  by  appropriate  affine  invariant  directions. 

Remark  4.11.  If  the  triangulation  covers  only  a  compact  domain  E.  then  usual  modifica¬ 
tions  of  basis  functions  corresponding  to  boundary  edges  or  vertices  (see  [22,  23])  lead  to  the 
desired  stable  local  bases. 

Remark  4.12.  In  this  section,  we  extended  to  the  setting  of  SLR-triangulations  the  bivari¬ 
ate  version  of  nodal  stable  local  basis  construction  of  [22,  23],  which  was  originally  designed 
for  regular  triangulations.  The  scheme  from  [27]  can  be  used  as  an  alternative  means  of 
constructing  stable  local  bases  for  Sk'r(Tm),  k  >  3 r  +  2,  in  dimension  cl  =  2.  Such  a  devel¬ 
opment  would  take  advantage  of  the  affine  invariance  of  the  Bernstein-Bezier  representation 
of  piecewise  polynomials.  We  elected  to  utilize  the  scheme  from  [22]  instead,  since  it  is 
available  for  any  number  of  variables  and  allows  an  effective  numerical  implementation  as 
shown  (for  r  =  1,2,  cl  =  2)  in  [23].  Also,  we  want  to  pay  heed  to  two  more  spline  basis 
constructions  (for  regular  triangulations  in  dimension  cl  =  2)  that  allow  the  same  kind  of 
extension  to  SLR  triangulations:  (a)  stable  local  bases  for  Sk'l(Tm, ),  k  >  5,  constructed 
in  [26]:  (b)  locally  stable  bases  on  nested  triangulations  (k  >  4 r  +  1)  [24],  Note  that  the 
stable  local  bases  for  superspline  subspaces  of  Sk,r(Tm )  [16,  17,  44,  57]  cannot  be  used  since 
these  spaces  are  not  nested  for  nested  triangulations,  while  the  earlier  local  spline  bases  for 
Sk,r(Tm)  [1,  8,  18,  35,  36,  48]  are  known  to  be  unstable  for  certain  triangulations. 

Remark  4.13.  It  is  easy  to  see  that,  in  the  case  r  =  0,  the  above  basis  reduces  to  the  classical 
Lagrange  finite  element  basis  for  Skfi(Tm),  k  >  1.  Since  82  disappears  from  (4.8)  when  r  =  0, 
Theorems  4.4  4.6  hold  for  locally  regular  triangulations,  i.e.,the  SLR  assumption  (2.3)  is  not 
needed  in  this  case.  (Note  that  S2  and  N0  completely  disappear  from  Theorem  4.4,  and  S2  is 
replaced  by  di  in  Theorems  4.5  4.6.)  For  r  =  0,  k  =  2,  we  get  the  Courant  elements,  and  the 
only  essential  difference  to  the  construction  from  [38]  is  that  we  rely  here  on  the  extensions 
of  linear  functionals  described  in  Remark  4.1  rather  than  on  the  explicit  quasi-interpolant 
for  continuous  piecewise  linear  functions  adopted  in  [38].  Both  approaches  obviously  lead  to 
the  same  B-spaces. 

5  Spline  bases  on  special  triangulations 

There  are  several  constructions  of  differentiable  spline  bases  fitting  into  our  scheme  that 
are  only  available  for  specific  multilevel  triangulations.  Since  these  triangulations  have  a 
special  structure  or  even  are  uniform,  the  corresponding  libraries  {Tr}  °f  bases  are  not 
as  rich  as  the  one  of  the  previous  section  associated  with  arbitrary  SLR-triangulations. 
Moreover,  the  necessity  to  maintain  the  structure  of  the  triangulation  highly  reduces  the 
variety  of  refinement  methods  that  can  be  used  (whereas,  e.g.,  local  refinement  by  bisection 
can  be  used  with  bases  on  arbitrary  triangulations.)  On  the  other  hand,  bases  on  special 
triangulations  usually  allow  a  smaller  degree  of  piecewise  polynomials  for  a  given  order  of 
differentiability  as  well  as  a  simpler  and  more  efficient  practical  implementation. 
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In  this  section,  we  review  some  known  constructions  of  this  type.  (Note  that  only  box 
splines  are  available  for  more  than  two  variables.) 

5.1  Box  splines 

As  usual,  we  consider  only  splines  of  two  variables.  Let  E  =  [£i  •  •  •  £n]  be  a  full  rank  2  x  n 
matrix  with  columns  C;  in  Z2  \  0.  The  box  spline  Me  :  R2  — >■  R2  associated  with  E  is  defined 
by  its  Fourier  transform 


mu-, = p 

i/=i 


i 


_  g  -i£vu 

iivU 


u  G  R2 , 


where  £vu  denotes  the  inner  product  of  the  two  vectors. 

We  now  review  the  basic  properties  of  box  splines  (see  [9]),  in  order  to  verify  the  require¬ 
ments  of  §  2.2.  It  is  well-known  that  Ms  has  a  compact  support, 

n 

supp  Me  =  |  ■  o  <  tv  <  1  j.  (5.1) 

v=\ 

The  box  spline  basis  functions  at  the  mth  level  are  defined  by 

¥Vj  =  MM1"'  •  -  j),  j  G  Z2. 

We  set 

I’m  =  {frj  :  i  6  Z  } ,  WJ  G  Z, 

and 

Sm.  —  ^  ^  ^  ■  &m,j  £  Rj-  i  ^ 

je  z2 

where  the  series  converges  everywhere  since  for  every  x  G  R2  and  m  G  Z  only  a  finite  number 
of  (j  G  Z2)  are  nonzero.  Clearly,  any  affine  change  of  variables  Q  :  R2  — >■  R2  gives 

rise  to  basis  functions  pm,j{Qx)  that  satisfy  the  conditions  of  §  2.2  if  and  only  if  the  do. 
Therefore,  we  do  not  distinguish  between  constructions  that  can  be  transformed  into  each 
other  by  such  a  method. 

Since 

Ms(2u) 

Me{u) 

Me  is  a  finite  linear  combination  of  M=( 2  •  —  j),  j  G  Z2,  which  implies  that 


n 

V=1 


1  +  e 


—i^vU 


Sm  C  <Srra+l,  G  Z. 


Let 

r(£)  :=  max{r  :  any  2  x  (n  —  r)  submatrix  of  E  has  rank  2}  —  1 

and 

A:(E)  :=  n  —  1. 
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The  elements  of  Sm  are  /’(E)  times  differentiable  piecewise  polynomials  of  degree  A:(S)  —  1 
with  respect  to  the  rectilinear  partition  Tff  of  R2  determined  by  the  straight  lines 

H„  +  2~mj,  jTZ2,  u  =  l  ,...,n, 


where 

Hv  :=  \tfv 

:  t  G  M} . 

Thus, 

Sm  C  Sk&- 

Moreover, 

hfc(=)  C  Sm, 

m  G  Z, 

and  n*(s)+1  (jL  Sm.  where 

A:(E)  =  r(S)  +  2. 


It  is  well-known  that  the  translates  of  a  box  spline  are  not  always  linearly  independent. 
In  fact,  <hm  is  a  basis  for  Sm  (m  €  Z)  if  and  only  if  the  matrix  E  is  unimodular ,  i.e.,each 
nonsingular  2x2  submatrix  of  E  has  determinant  ±1.  This  condition  implies  substantial 
restrictions  on  E.  Namely,  up  to  an  affine  change  of  variables,  E  must  have  the  form 


where  e3  =  [J] ,  e2  =  [°] ,  e3  =  []] ,  n3,  n2  >  1,  n3  >  0,  and  +  n2  +  n3  =  n.  It  is  easy  to 
see  that 

r(E)  =  n  —  max{ni,  n2,  n3}  —  2 

and  that  is  either  a  tensor  product  mesh  if  n3  =  0  or  a  three- directional  mesh  Tin'1  defined 
by  the  straight  lines  x\  =  2 ~mj,  x2  =  2 ~mj,  x\  —  x2  =  2 ~mj  ( j  6  Z2)  in  M2  if  n3  >  1.  Since 
only  the  latter  case  leads  to  a  multilevel  triangulation,  we  assume  that  n3  >  1. 

It  remains  to  verify  (2.10)  (2.12).  By  (5.1),  the  support  of  Ms  is  the  hexagon  with 
vertices  (0,0),  (ni,0),  (0,  n2),  (n3  +  n3,n3),  (n3,  n2  +  n3),  (n3  +  n3,n2  +  n3),  which  implies 
(2.10)  with  A  <  1_7t./2J  .  Obviously,  (2.11)  is  valid  with  M3  =  ||Ms||Loo.  Finally,  it  is  easy 
to  show  (2.12)  by  using  the  constructions  of  dual  functionals  A  j  :  <So  — >■  M  (j  G  Z2),  with 
-M^cufc)  =  given,  e.g.,  in  [19,  37,  41]. 

Let  us  mention  the  following  two  cases  that  are  perhaps  most  relevant  in  applications: 

(a)  m  =  n2  =  2,  n3  =  1,  Sm  C  S4'1  {Tm] ) t  A:  =  3, 

(b)  m  =  n2  =  n3  =  2,  5m  C  <S5,2(Tm '),  k  =  4. 

5.2  Other  spline  bases  on  uniform  triangulations 

There  are  some  other  spline  basis  constructions  for  the  three-directional  mesh  Tm  \  see,  e.g., 
[15,  56].  However,  to  our  knowledge,  none  of  them  simultaneously  satisfies  the  requirements 
of  nestedness  of  the  spaces,  stability,  and  locality  of  the  basis  functions.  The  situation  is 
better  for  the  four- directional  mesh  Tn?  obtained  from  Tm'1  by  adding  the  straight  lines 
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X\  +  x2  =  2 ~mj  ( j  G  Z2).  Since  Tm'1  is  a  special  case  of  a  so-called  FVS-triangulation  (see 
Section  5.3),  finite  element  bases  for  S4'1(Tm'>)  are  available  and  satisfy  the  conditions  of 
Section  2.2.  Some  recent  alternative  constructions  of  stable  local  bases  for  <S4J  (Tm  ) )  can  be 
found  in  [13,  28,  42,  49].  Moreover,  a  stable  local  basis  for  S7-2(Tn?)  is  also  constructed  in 
[28] .  Finally,  we  want  to  mention  the  stable  local  basis  from  [33]  for  C1  quadratic  splines  with 
respect  to  a  sequence  of  triangulation  levels  that  can  be  called  the  six-directional  meshes. 

5.3  Refinable  composite  finite  elements 

Multilevel  and  hierarchical  bases  play  an  important  role  in  the  modern  theory  and  practice 
of  numerical  methods  for  PDEs,  see,  e.g.,  [51].  Classical  smooth  finite  elements  [14]  give  rise 
to  stable  local  spline  bases  on  triangulations  satisfying  the  minimal  angle  condition.  (Note 
that  it  should  be  possible  to  replace  this  condition  of  regularity  with  SLR.)  However,  there 
are  difficulties  in  using  them  to  build  nested  spline  spaces  on  multilevel  triangulations,  see 
[11,  20].  Although  the  “polynomial”  finite  elements  (e.g.,  Argvris  element)  are  available  for 
arbitrary  triangulations,  they  lead  to  superspline  spaces  [57]  that  lack  nestedness  for  nested 
triangulations  ( levels  in  the  terminology  of  our  §  2).  In  contrast  to  them,  “composite” 
finite  elements  require  a  special  structure  of  the  levels  Tm ,  e.g.,  a  Clough-Tocher  or  Powell- 
Sabin  split,  which  is  not  always  compatible  with  nested  refinements  with  other  desirable 
properties  like  boundedness  of  the  valence  of  the  vertices.  In  fact,  we  are  aware  of  only  two 
cases  when  composite  finite  elements  are  refinable,  i.e. ,  provide  stable  local  bases  for  certain 
multilevel  triangulations.  First,  this  is  true  for  the  triangulations  obtained  by  the  Powell- 
Sabin  12-split,  see  [50]  for  the  relevant  construction  of  stable  local  bases  for  C1  quadratics 
and  cubics.  The  other  case  is  that  of  FVS-triangulations  obtained  from  arbitrary  srictlv 
convex  quadrangulations  by  adding  two  diagonals  of  each  quadrilateral,  see,  e.g.,  [20,  43]. 
Here,  a  well-known  composite  finite  element  due  to  Fraeijs  de  Veubeke  and  Sander  gives 
rise  to  a  stable  local  basis  for  C 1  cubics,  while  for  higher  orders  of  differentiability  only 
non-nested  superspline-type  constructions  are  known  [40,  45,  46]. 


6  Appendix 

Proof  of  Theorem  2.9.  Denote  briefly  N  :=  (^ 
Case  1 :  0  <  p  <  1.  Since  t  <  p  <  1,  we  have 


ic^e(-)i 

0e@ 


< 


(Elfins) 

oe@ 


0e& 

i/p 


ll<w«l  pv'> 

<  (Eiiwii;) 1  ■ 

0ee 


Case  2:  p  =  oo.  Since  r  =  1/a  <  1,  then  (2.42)  is  obvious. 

Case  3:  1  <  p  <  oo.  We  need  the  following  lemma: 

Lemma  6.1.  Let  g  :=  \C0LP0 \,  where  #A4  <  oo  and  He^^Hp  <  L  for  9  G  A4.  Then 

llsllp  < 


with  e  independent  of  M.  and  {ce)eeM- 


44 


Proof.  Using  the  properties  of  T,  we  have  (recall  that  supp^<?  C  Eg  :=  star  e(vg)  and 
Halloo  ~  \Ee\-^p\\(p0\\p  by  (2.14)) 


IIA<  E  ||c0V?e[|oo  •  Ifie(-) 


<  cL\ 


oeM 


E 

OeM 


E, 


I  -i/p 


11 /:«(') 


We  define  E  :=  \J0eM  Eg  and  E(x)  :  =  min{|E0|  :  0  E  M  and  Eg  3  x }  for  x  €  E.  By  the 
properties  of  the  LR-triangulations,  it  follows  that 


5]  \Eg\-^  ■  lEe(x)  <  cE(x)~l/nE{x),  X  e  K2. 
OeM 


On  the  other  hand, 


E(x)  *=  \Eo\  1  <  \E»\  llEetx)- 

OeM ,  Eg3x  z ' 


OeM 


Therefore, 


Mp  <  cL[|i?(-)-1/p|UP  =  cL(^E(.r)-Vte)  ’ 

<  CL(  £  I  Eel-1  f  1  Eg(x)  dx)1/P  =  cL(#M)1/E 
oeM 


We  define 


Tp  :=  {0  :  2  "N  <  \\co<pB\\p  <  2^+1N}, 


where  N  :=  {J^oee  \\coVe\\l)llr ■  Then 


U  Tv  =  {e  :  WweWp  >  2"^} 

U</J: 


and  hence 

=  #(  U  <  2"r- 

v</i  ^</J, 

We  set  F/(  :=  YhoeXg  Using  Lemma  6.1  and  (6.1),  we  obtain 

OO  OO  OO 

||Ei‘w»(-)i  ^  II  Ew  <Eii^ii-^<!E2"''n(#g.)1/!' 


(6.1) 


0e& 


fj— o  X  p=0  p=0 

OO  OO 

<  cN  ^  2-p(1-r/p)  <  cN  ^  2~pra  <  cN. 

l-i— 0  i-i—O 


This  completes  the  proof  of  Theorem  2.9. 


Proof  of  Theorem  3.1  (Case  0  <  p  <  oo).  Suppose  /  €  £?“($),  where  a  >  0, 
1/r  =  ck  +  1/p,  0  <  p  <  oo.  By  (2.40),  /  can  be  represented  in  the  form  /  =  ^20ee  bgp)g  with 
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the  series  converging  absolutely  a.e.  in  R2  and  in  Lp.  We  denote  briefly  N(f)  :=  :  = 

(E0g©  \\bo^e\\Tp)1/r  ~  I  ./'I  /»•■:•  Cl'S- 

Suppose  that  (bg  )J  is  a  rearrangement  of  the  sequence  {bgipg)geQ  such  that  > 

1 1 6(92  W02  Up  >  ■  ■  ■■  Set  sn  :=  Y,'-  j  A,  ro,  •  s„  G  £„($). 

Case  /:  0  <  p  <  1.  To  estimate  [|/  —  sn||p  we  shall  use  the  following  simple  inequality 
[38] :  If  x\  >  X2  >  •  •  •  >  0  and  0  <  r  <  p,  then 

(E*5)1/'s«,/-1/'(f:*;)1/T.  (6-2) 

./  M  •  I  ./  I 

We  use  Theorem  2.9  and  apply  (6.2)  with  Xj  :=  \  bg  | to  obtain 

0°  OO  -j Jp 

\\f-sn\\p  <  |  ^  ^  WbojWjW l) 

j=n-\- 1  j=n+ 1 

OO  H  / 

<  «1/^1/T(£llM«iK)  =n-“JV(/), 

J=1 

which  proves  Theorem  3.1  in  Case  1. 

Case  2:  1  <  p  <  oo.  We  proceed  quite  similarly  as  in  the  proof  of  Theorem  2.9.  We  set 
T,  :=  {6  :  2  ",Y(/!  <  ||6^||p  <  2 ~»+1N(f)}  and  Fp  :=  |6,<pfl|. 

Fix  m  >  1  and  set  M  :=  |~2TOT~|.  As  in  the  proof  of  Theorem  2.9  (see  (6.1)),  < 

<  2mr  <  M.  Using  Lemma  6.1,  we  obtain 

OO  OO 

||/  —  sm|Ip  A  ^  ^  "  ll-f/llp 

p 

/x=m+l  /i=m+l 

oo  oo 

<  c  J]  2 -C\(f)(#F,)^<cN(f)  2-»(1-T/p) 

H=m+ 1  n=m+ 1 

<  cW(/)2-m(1-T/p)  <  r.l/  1  •  1  CV(/!  =  cM~aN(f). 

This  estimate  readily  implies  (3.1).  ■ 

Proof  of  Theorem  3.2.  Step  1.  With  this  step  we  lay  some  groundwork  that  is  needed  for 
the  proof  of  the  Bernstein  inequality.  Let  T  be  an  arbitrary  LR-triangulation  and  suppose 
A  is  a  finite  subset  of  T.  The  set  A  generates  a  certain  tree  structure  that  we  want  to  bring 
up  in  what  follows. 

We  say  that  A  e  T  is  a  branching  triangle  if  at  least  two  children  of  A  have  descendants 
in  A.  Let  A  denote  the  extension  of  A  obtained  by  adding  all  branching  triangles  and  all 
children  of  branching  triangles  if  they  are  not  already  in  A.  By  considering  the  tree  of  the 
ancestors  of  all  triangles  in  A,  it  is  not  difficult  to  see  that  the  total  number  of  branching 
triangles  does  not  exceed  #A  —  1.  Since  the  number  of  children  of  a  triangle  is  bounded  by 
M0,  we  conclude  that  #A  <  c# A. 

Furthermore,  for  a  later  use  in  Step  3,  we  call  A  e  T  \  A  a  chain  triangle  if  at  least 
one  of  its  descendants  belongs  to  A.  The  set  of  all  chain  triangles  will  be  denoted  by  F.  By 
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construction,  for  each  A  6  F  there  is  a  unique  largest  triangle  A  G  A  contained  in  A.  We 
set  /\'a  :=  A  \  A  and  //a  :=  m  —  m,  where  A  G  7^  and  A  e  7^,,  We  denote  by  F  the 
set  of  all  A  e  T  for  which  there  is  a  A'  €  A  containing  A.  It  is  easy  to  see  that  F  is  the 
disjoint  union  of  finite  chains ,  i.e. ,  sets  A  of  the  form  A  =  {Ax, . . . ,  A„}  C  T  (v  >  1),  where 
A"  D  Ax  D  ■  ■  ■  D  A„  D  A'x  for  some  A'a,  A"  G  A,  and  Ax  is  a  child  of  A",  A?  is  a  child  of 
Aj_x,  v  =  2, . . , ,  v,  and  A'x  is  a  child  of  A„.  Similarly,  T  \  T  is  the  disjoint  union  of  infinite 
chains  A  =  {. . . ,  A_2,  A_x}  C  T,  where  •  •  •  D  A_2  D  A_x  D  A'x  for  some  Aa  G  A,  and  A j 
is  a  child  of  Aj_x,  v  =  —1,  —2, . . .,  and  A'x  is  a  child  of  A_x.  We  let  C  and  C°°  denote  the 
sets  of  all  finite,  respectively,  infinite  chains  in  T.  Clearly,  ffC  <  fiA  and  ffC°°  <  #A. 

Step  2.  For  the  proof  of  the  theorem  in  the  case  0  <  p  <  oo,  we  need  the  following 
lemma. 


Lemma  6.2.  Suppose  s  =  ^ag4^a  '  where  Pa  £  (k  >  1).  Ac  T  with  T  an 
LR-triangulation,  and  fiA  <  oo.  Then 

( J2  iAr”isiii#^)1/T  <  c(#a)°iis|Ip 

aga 


with  c  independent  of  s  and  A. 

Proof.  We  adopt  all  necessary  notation  from  Step  1  above  with  A  from  the  hypotheses  of 
the  lemma.  Since  #A  <  c# A  and  s  =  Yl/s..  \  '  A  -  where  Pa  =  0  for  A  €  A  \  A,  we  may 

assume  without  loss  of  generality  that  A  =  A,  i.e.,  the  branching  triangles  and  their  children 
are  contained  in  A. 

Let  Ax, ... ,  Am  be  all  non-branching  triangles  in  A.  It  is  not  difficult  to  see  that  for  each 
of  them  there  are  only  two  possibilities:  either  A.,  does  not  contain  any  other  A  6  A  (in 
which  case  we  call  A,;  a  fi,nal  triangle )  or  there  is  a  unique  largest  triangle  A,;  €  A  strictly 
contained  in  A*.  We  define  the  rings  Kj  :=  A*  \  Aj,  i  =  1, . . . ,  m,  where  A*  :=  0  for  a 


final  triangle  Aj 

P  c  nA:.  i  =  i, . 

A  G  A, 


Obviously,  I(,  have  pairwise  disjoint  interiors,  and  s |  /x,  =  P, |  k,  for  some 
. ,  m.  Since  all  children  of  branching  triangles  are  in  A,  we  have  for  each 


A  = 


Therefore, 


EGI 

aga 


lir(A) 


m  m 

U  Ki  and  s|A  =  tKi  •  P- 

l=L 

AiCA 

AjCA 

m 

.  =  E 

|A|—  E  WILw 

AgA 

i=l 

AjCA 

m 

=  Ei 

MILm  E  lAl_“T 

2=1 

A£A,  Ad  A; 

m 

=  Ei 

l*IIL(,fi)|A.r“T  E 

2=1 

AeA,  ADAj 

m 

(|A?|/|A|r 


< 


E 

2=1 


/.q/'i)  A; 
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where  we  once  switched  the  order  of  summation  and  used  (2.9).  Since  |A,;|  <  p|A,;|,  we  have 
by  Lemma  2.2, 


i\\LT(Ki) 


\KA1/t~1/p\\R; 


i\\Lp(Ki) 


I  A,; 


i\\Lp{Ki), 


which  implies  ||  s  1 1  I  A*  |  ar  ~  ||s||l  (^p  *  =  ,m.  Now  by  Holder’s  inequality, 


m 


i— 1 


m  , 

(EiisHA,))T'’mW'<(#Arik 


and  the  proof  is  complete.  ■ 


Step  3.  Let  s  G  En(<f>)  and  suppose  that  s  =:  w'i ere  -M.  C  0(T)  and 

#A4  <  n.  Let  A  be  the  set  of  all  triangles  A  G  T  which  are  involved  in  all  Eg  :=  supp  pg, 
6  G  M.  Then  s  =  SM  w^ere  -sa  =:  1a  •  -Pa,  Pa  €  II*.  Evidently,  by  (2.7),  # A  < 

A  (A',,./)  /I  M  <  <)i. 

We  first  extend  A  to  A  as  in  Step  1  above  and  introduce  some  auxiliary  sets  of  triangles 
needed  for  the  forthcoming  arguments.  We  set 

A*n  :=  {A  G  Tm  ■  11/x  D  A7  for  some  A'  G  A  n  Tm}, 

A**  :=  e  7 m  ■  I^a  ^  H 7  for  some  A7  G  A  fl  Tm},  m  G  Z, 

and  also 

A*  :=  U  A  A”  :=  [J  A” 

mG  Z  me  z 

Note  that  A,  A7  G  7^.  and  A7  C  imply  A  C  and  hence 

AW  {A  G  Tm  :  A  C  for  some  A7  G  A  fl  7^}. 

Therefore,  by  (2.8),  #A*n  <  c**(iV0,  t)  #(A  fl  7^),  and  it  follows  that  #A*  <  cn.  Similarly, 
#A**  <  c**(iV0,  2(y)(#A)  <  cn.  It  is  clear  that  A  C  A*  C  A**. 

We  now  proceed  to  estimate  :=  J^Aer  I  A|_qtSa(s)^.  Let 

sm  ■=  ^2  c° V0'  m  G  Z- 

\i<m  oeMn&fi 

Then  sm  G  Sm  and  hence  §A(s)T  =  SA(s  —  sm)T  if  A  G  Tm.  For  each  A  G  T,  we  shall  use 
one  of  the  following  two  obvious  bounds  for  SA (s)T: 


IIsIIlt(o<a), 

(6.3) 

IIs  —  •c’ih  1  / .  ;o\}' 

A  G  Tm  . 

(6.4) 

Namely,  (6.3)  will  be  applied  to  the  triangles  A  in  the  set  A*  C  T  defined  above,  while  (6.4) 
will  be  used  for  all  remaining  triangles  in  T. 

For  the  next  estimates,  we  shall  consider  separately  the  cases  0  <  p  <  oo  and  p  =  oo. 
Case  1 :  0  <  p  <  oo.  We  consider  two  possibilities  for  each  A  G  T:  A  G  A*  or  A  G  T\  A*. 
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(a)  If  A  G  A;,,  then  for  each  A'  G  Tm  such  that  A'  C  QeA,  we  have  A'  G  A”  and,  in  view 
of  (2.2),  | A7 1  <  c|A|.  Hence,  by  (6.3), 

5]  |A|— S A«  <  V  |A|—  5]  ll«IIL(A') 


<  <E 

E  iA'r°Tt 

IsIIlt(a') 

A'GA-,A'cOi 

II 

M 

E  lA'l"G 

lsilI,T(A') 

a'ga  ” 

AeA*m,Q(ADA' 

£  *  E 

Ia'|-“tI|S||At 

where  we  have  switched  the  order  of  summation  and  taken  into  account  the  fact  that  #{A  G 
:  Q£a  D  A'}  =  #{A  G  A^  :  A  C  12^,}  <  c**(IY0,  £),  by  (2.8).  It  follows  that 


5]  |A|— sA(S);  <  c  5]  |Ar“T||s||y(A) 

AeA*  AgA** 

<  c(#A**)"T||s||p  <  cnaT||s||p,  (6.5) 

where  we  applied  Lemma  6.2  to  s  with  A  replaced  by  A**,  which  is  obviously  legitimate  since 
A**  d  A. 

(b)  Now  suppose  A  G  Tm\Kn ■  Then  QfA  =  Uj=i  A j  for  some  Aj  G  Tm  \A,  j  =  1 . iia. 

with  nA<  A*  <  3Nq£_1  (see  (2.8)).  We  have,  using  (6.4), 

«A 

Sa(s);  =  S  a(s  -  sm)rT  <  ^  |  s  S„,\  (6.6) 

3= 1 


If  Aj  j.  T,  then  it  has  no  descendants  in  A,  and  hence  s\A.  =  sm\Aj,  and 

II  ^  1 1  £rff  ( A  j1)  6 ,  Aj  <f  P. 


(6.7) 


Suppose  Aj  G  P,  i.e. ,  it  is  a  chain  triangle.  Let  A  j  be  the  unique  largest  triangle  of  A 
contained  in  Aj,  and  let  KAj  =  Aj\Aj  and  [iAj  =  m  —  fh  be  defined  as  in  Step  1.  It  is  clear 
that  in  this  case  s\Ka.  =  sm \Ka  =  1Ka.  ■  PAj  and  sm\Aj  =  lAj  ■  PAj  for  some  PAj  G  Uk. 
Therefore, 


S  -  Sr 


l.T  ( Aj ) 


=  S  —  Sr 


(A)  - 


<  ells 


Aj) 


c\\Pt 


A,- 


r 

Lr(Aj) 


If  Aj  G  T  \  T,  then  clearly  =  0,  and  we  have 


\s  - 


I  -^r(Aj) 


=  \\s\ 


Lr(Aj)  ’ 


Aj  e  r  \  f . 


(6.8) 


Assume  that  Aj  G  P.  By  Lemma  2.2, 


11^.  IIlt(aj) 


< 

< 


\A\\PaMaA)  -  Cl^illl-PA,-||l00(A'Aj) 

C|Aj||/iAj|-^||^j||Ip(sv<<i|Aj||Aj|"-1||i!||Ii(jrii.p? 


49 


By  (2.1),  |Aj|  <  |  A j | ,  and  we  arrive  at  the  inequality 

IIs  -  Sm\\TLT(Aj)  <  441^)  +  A/,/|.s'|';(,KA(,. 

From  (6.6)  (6.9)  and  (2.2),  we  obtain 


A,  e  r. 


E  iai““ts*m; 

4«t\a* 


=  E  E  iai_“tsam; 

m&  AeTm\A*m 

<  cEiai-'NII«+«EowMm^ 


Aer 

Si  +  S2. 


Aer 


Trivially, 


Ilt(A)  <  E  llS'l|ir(A'),  Aer. 
A'eA,  A'C A 


Switching  the  order  of  summation,  we  find 


< 

A 

I  s| 

L(A') 

E 

|A|-“T 

A'eA 

Aer,  Ar  A' 

< 

<E 

I  s| 

IL(A') 

b- 

3 

i 

E  (iA'i/iAir 

A'eA 

Aer,  Ada' 

< 

cE 

|A‘ 

'|—  ar  ||  ( 

’i  L(A')  < 

c(#A)“T||s||;, 

A'eA 

(6.9) 


(6.10) 


where  we  also  used  (2.9)  and  applied  Lemma  6.2  to  s  with  A  replaced  by  A. 

To  estimate  S2  we  shall  use  the  representation  of  F  as  a  disjoint  union  of  chains:  T  = 
|JAe£  A.  Let  A  e  C  and  suppose  A  =  {Ax, . . . ,  A„},  where  A"  )  A)  )  •  •  •  )  A„  D  Aa  with 
A'a,  A"  e  A.  Then  //a,  >  v  —  i  +  1.  Therefore, 


E 

AeA 


Pij'a\\s\\lp(ka)  E  HsIIl1)(a"\a'a) 


<  llsl 


E  py~1+l  < 


c  s 


Lp(K\)  ’ 


i=  i 


where  Kx  :=  A"  \  Aa.  It  is  easy  to  see  that  the  sets  A  A  •  A  G  C.  have  pairwise  disjoint 
interiors.  Summing  over  all  A  G  £,  we  obtain  by  Holder’s  inequality 


A  <  cE  IMIG,™  <  AE  llsllh(M.)T/’’(#£)‘-T/*'  <  C(#A)-||S 


A ec  A e£ 

From  this  estimate  and  (6.10),  we  find 


E  |a|-“tsa(s);  <  c(#A) 

Aer\A* 


“™s|£  <  cnaT\\s\\Tp. 


Combining  this  with  (6.5)  gives  ||s||gQ^  <  cn"T||s||p,  i.e. ,  (3.2)  holds. 
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Case  2:  p  =  oo.  The  proof  in  this  case  is  easier.  We  consider  as  before  two  possibilities 
for  each  A  G  T:  A  G  A*  or  A  G  T  \  A*. 

(a)  For  A  G  A*,  we  obtain  by  (2.2) 

|A|_1Sa(s);  <  |A|-1||,|rMQi)  <  | A|_1|fi^|||s||J0  <  c||,||^. 

Therefore, 

E  |A|_1§a(S);  <  c||s|iJ0(#A*)  <  HI*-  (6.11) 

AeA* 

(b)  Let  A  G  Tm  \  A *m.  Then  =:  Uj=i  f°r  some  A j  G  Tm  \  A,  j  =  1, . . . ,  //A;.  with 
nA  <  c**  <  3Aq£_1  (see  (2.8)).  We  have  (see  (6.4)) 

nA 

SA(s);  =  §A(s  -  STO)r  <  Y  IIs  -  SmWlgQ,)- 

3= 1 

As  in  Case  1 ,  if  Aj  ^  T,  then  ||.s  —  £to||lt(Aj)  =  0,  and  if  A  j  G  T,  then  s\Ka  =  sm\KA  = 
1  ka.  ■  Paj  and  sm |Aj.  =  lAj  •  PAj  for  some  PA.  G  II*.  Therefore, 

IIs  “  Srn\\LT(Aj)  =  IIs  —  /  .  i 

<  +  IIPaJI^)  <  -  a,|.s-k. 

where  we  used  the  inequalities  I  Pa,  \  ,  N  0j}  <  I  Pa  ,  \  usa,i  <  c||PaJ£<30(a-a  ,  <  c|  *|  x  (see 
Lemma  2.2).  From  the  above,  we  infer  by  (2.2) 

|A|-‘sa(S);  <  c||Sfc  E  IGI/IGI. 

Aj-ernTm 

and  hence,  using  (2.2)  and  the  fact  that  each  A'  G  T  n  7^,  can  belong  to  <  c**  sets  we 
obtain 

E  |Ar'Sa(S);  <  c\\s\\l  E  |A|/|A| 

AeTm\Arn  Aero?;,, 

<  «jisii^o  Y  p‘JA- 

Aer  nTm 

Summing  over  m  G  Z,  we  find 

E  I AI-'SaMt  <  E  !f“  £  cll*ll«  (#c  +  £  “MIL- 

AeT\A*  Aer 

We  couple  this  with  (6.11)  to  obtain  ||s||gQ^  <  cnHsH^,  which  is  (3.2).  ■ 


51 


References 


[1]  P.  Alfeld,  B.  Piper,  and  L.  L.  Schumaker,  Minimally  supported  bases  for  spaces  of 
bivariate  piecewise  polynomials  of  smoothness  r  and  degree  d  >  4 r  +  1,  Comput.  Aided 
Geom.  Design,  4  (1987),  pp.  105  123. 

[2]  I.  Babuska  and  A.  K.  Aziz,  On  the  angle  condition  in  the  finite  element  method ,  SIAM 
J.  Numer.  Anal.,  13  (1976),  pp.  214  226. 

[3]  J.  Bergh  and  J.  Lofstrom,  Interpolation  spaces.  An  introduction ,  Grundlehren  der  Math- 
ematischen  Wissenschaften,  No.  223.  Springer- Verlag,  Berlin-New  York,  1976. 

[4]  C.  Bennett  and  R.  Sharpley,  Interpolation  of  operators ,  Pure  and  Applied  Mathematics, 
129.  Academic  Press,  Inc.,  Boston,  MA,  1988. 

[5]  P.  Binev  and  R.  DeVore,  Fast  computation  in  adaptive  tree  approximation ,  preprint, 
(http:  / /www.math.sc.edu/^imip/02.html) 

[6]  P.  Binev,  W.  Dahmen,  and  R.  DeVore,  Adaptive  finite  element  methods  with  convergence 
rates ,  IGPM  Report  No.  219,  RWTH  Aachen,  June  2002. 

(http:  / /elc2. igpm.rwth-aachen.de/ ^dahmen/) 

[7]  P.  Binev,  W.  Dahmen,  R.  DeVore,  and  P.  Petrushev,  Approximation  classes  for  adaptive 
methods ,  Serdica  Math.  J.,  28  (2002),  pp.  391  416. 

[81  C.  de  Boor  and  Iv.  Hollig,  Approximation  power  of  smooth  bivariate  pp  functions ,  Math. 
Z.,  197  (1988),  pp.  343  363. 

[9]  C.  de  Boor,  K.  Hollig,  and  S.  D.  Riemenschneider,  Box  Splines,  Springer- Verlag,  New 
York,  1993. 

[10]  C.  de  Boor  and  R.  Q.  Jia,  A  sharp  upper  bound  on  the  approximation  order  of  smooth 
bivariate  pp  functions ,  J.  Approx.  Theory,  72  (1993),  pp.  24  33. 

[11]  J.  Bramble  and  X.  Zhang,  Multigrid  methods  for  the  biharmonic  problem  discretized  by 
conforming  C1  finite  elements  on  nonnested  meshes ,  Numer.  Funct.  Anal.  Optim.,  1 
(1995),  pp.  835  846. 

[12]  S.  Brenner  and  L.  R.  Scott,  The  mathematical  theory  of  finite  element  methods ,  Texts 
in  Applied  Mathematics,  15,  Springer- Verlag,  New  York,  1994. 

[13]  M.  D.  Buhmann,  O.  Davydov,  and  T.  N.  T.  Goodman,  Cubic  spline  prewavelets  on  the 
four- directional  mesh ,  Found.  Comp.  Math.,  3  (2003),  pp.  113  133. 

[14]  P.  G.  Ciarlet,  The  Finite  Element  Method  for  Elliptic  Problems ,  Nortli-Holland,  Nether¬ 
lands,  1978. 

[15]  C.  Iv.  Chui,  Multivariate  Splines,  SIAM,  Philadelphia,  1988. 


52 


[16]  C.  Iv.  Chui,  D.  Hong,  and  R.-Q.  Jia,  Stability  of  optimal  order  approximation  by  bivariate 
splines  over  arbitrary  triangulations ,  Trans.  Amer.  Math.  Soc.,  347  (1995),  pp.  3301 
3318. 

[17]  C.  K.  Chui  and  M.-.J.  Lai,  Multivariate  vertex  splines  and  finite  elements ,  J.  Approx. 
Theory,  60  (1990),  pp.  245  343. 

[18]  C.  K.  Chui  and  M.-.J.  Lai,  On  bivariate  super  vertex  splines ,  Constr.  Approx.,  6  (1990), 
399  419. 

[19]  W.  Dahmen  and  C.  A.  Micchelli,  On  the  local  linear  independence  of  translates  of  a  box 
spline ,  Studia  Math.,  82  (1985),  pp.  243  262. 

[20]  W.  Dahmen,  P.  Oswald,  and  X.-Q.  Shi,  Cl -hierarchical  bases ,  J.  Comput.  Appl.  Math., 
51  (1994),  pp.  37  56. 

[21]  W.  Dahmen  and  R.  Stevenson,  Element-by- element  construction  of  wavelets  satisfying 
stability  and  moment  conditions ,  SIAM  .J.  Numer.  Anal.,  37  (1999),  pp.  319  352. 

[22]  O.  Davydov,  Stable  local  bases  for  multivariate  spline  spaces,  J.  Approx.  Theory,  111 
(2001),  pp.  267  297. 

[23]  O.  Davydov,  On  the  computation  of  stable  local  bases  for  bivariate  polynomial  splines , 
in  “Trends  in  Approximation  Theory”  (Iv.  Ivopotun,  T.  Lyclie,  and  M.  Neamtu,  Eds.), 
Vanderbilt  University  Press,  2001,  pp.  85  94. 

[24]  O.  Davydov,  Locally  stable  spline  bases  on  nested  triangulations ,  in  “Approxima¬ 
tion  Theory  X:  Wavelets,  Splines,  and  Applications”  (C.K.Chui,  L.L.Schumaker,  and 
J. Stockier,  Eds.),  Vanderbilt  University  Press,  2002,  pp.  231  240. 

[25]  O.  Davydov,  G.  Niirnberger,  and  F.  Zeilfelder,  Bivariate  spline  interpolation  with  opti¬ 
mal  approximation  order ,  Constr.  Approx.,  17  (2001),  pp.  181  208. 

[26]  O.  Davydov  and  L.  L.  Schumaker,  Stable  local  nodal  bases  for  C 1  bivariate  polynomial 
splines ,  in  “Curve  and  Surface  Fitting:  Saint-Malo  1999”  (A.  Cohen,  C.  Rabut,  and  L. 
L.  Schumaker,  Eds.),  Vanderbilt  University  Press,  2000,  pp.  171  180. 

[27]  O.  Davydov  and  L.  L.  Schumaker,  On  stable  local  bases  for  bivariate  polynomial  spline 
spaces ,  Constr.  Approx.,  18  (2002),  pp.  87  116. 

[28]  O.  Davydov  and  F.  Zeilfelder,  Scattered  data  fitting  by  direct  extension  of  local  polyno¬ 
mials  with  bivariate  splines ,  Advances  in  Comp.  Math.,  to  appear. 

( http:  // www.uni-giessen.de  /  ~  gcn5  /  davydov  / ) . 

[29]  R.  DeVore,  B.  .Jawerth,  and  V.  Popov,  Compression  of  ivavelet  decompositions ,  Amer. 
J.  Math.,  114  (1992),  pp.  737  785. 

[30]  R.  DeVore,  P.  Petrushev,  and  X.  Yu,  Nonlinear  wavelet  approximation  in  the  space 
C(Rd ),  in  “Progress  in  Approximation  Theory”  (A.  A.  Gonchar,  E.  B.  Saif,  Eds.),  New 
York,  Springer- Verlag,  1992,  pp.  261  283. 


53 


[31]  R.  DeVore  and  V.  Popov,  Interpolation  spaces  and  nonlinear  approximation,  in  “Func¬ 
tion  spaces  and  applications  (Lund,  1986)”,  Lecture  Notes  in  Math.,  1302,  Springer, 
Berlin,  1988,  pp.  191  205. 

[32]  R.  DeVore  and  V.  Popov,  Interpolation  of  Besov  spaces ,  Trans.  Amer.  Math.  Soc.,  305 
(1988),  pp.  297  314. 

[33]  M.  Daehlen,  T.  Lvclie,  K.  Morken.  R.  Schneider,  and  H.-P.  Seidel,  Multiresolution  anal¬ 
ysis  over  trianqles,  based  on  quadratic  Hermite  interpolation ,  .J.  Comput.  Appl.  Math., 
119  (2000),  pp.  97  114. 

[34]  M.  Floater  and  E.  Quak,  Linear  independence  and  stability  of  piecewise  linear  pre¬ 
wavelets  on  arbitrary  triangulations ,  SIAM  .J.  Numer.  Anal.,  38  (2000),  pp.  58  79. 

[35]  D.  Hong,  Spaces  of  bivariate  spline  functions  over  triangulation.  Approx.  Theory  Appl., 
7  (1991),  pp.  56  75. 

[36]  A.  Ivh.  Ibrahim  and  L.  L.  Schumaker,  Super  spline  spaces  of  smoothness  r  and  degree 
d  >  3 r  +  2,  Constr.  Approx.,  7  (1991),  pp.  401  423. 

[37]  W.  Jianzhong,  On  dual  basis  of  bivariate  box  spline ,  Approx.  Theory  Appl.,  3  (1987), 
pp.  153  163. 

[38]  B.  Ivaraivanov  and  P.  Petrushev,  Nonlinear  piecewise  polynomial  approximation  beyond 
Besov  spaces ,  ACHA,  to  appear,  (http://www.math. sc. edu/^pencho) 

[39]  B.  Ivaraivanov,  P.  Petrushev,  and  R.  Sharpley,  Algorithms  for  nonlinear  piecewise  poly¬ 
nomial  approximation:  Theoretical  aspects,  Trans.  Amer.  Math.  Soc.,  355  (2003),  pp. 
2585  2631. 

[40]  M.  Laghchim-Lahlou  and  P.  Sablonniere,  Quadrilateral  finite  elements  of  FVS  type  and 
class  Cp,  Numer.  Math.,  70  (1995),  pp.  229  243. 

[41]  M.-J.  Lai,  A  remark  on  translates  of  a  box  spline.  Approx.  Theory  Appl.,  5  (1989),  pp. 
97  104. 

[42]  M.-J.  Lai,  Approximation  order  from  bivariate  C1- cubics  on  a  four- directional  mesh  is 
full ,  Comput.  Aided  Geom.  Design,  11  (1994),  pp.  215  223. 

[43]  M.-J.  Lai,  Bivariate  spline  spaces  on  FVS-triangulations ,  in  “Approximation  Theory 
VIII,  Vol.  1:  Approximation  and  Interpolation”  (C.  Iv.  Chui  and  L.  L.  Schumaker, 
Eds.),  World  Scientific,  Singapore,  1995,  pp.  309  316. 

[44]  M.-J.  Lai  and  L.  L.  Schumaker,  On  the  approximation  power  of  bivariate  splines ,  Ad¬ 
vances  in  Comp.  Math.,  9  (1998),  pp.  251  279. 

[45]  M.-J.  Lai  and  L.  L.  Schumaker,  On  the  approximation  power  of  splines  on  triangulated 
quadrangulations ,  SIAM  J.  Numer.  Anal.,  36(1)  (1999),  pp.  143  159. 


54 


[46]  M.-J.  Lai  and  L.  L.  Schumaker,  Quadrilateral  macro- elements ,  SIAM  J.  Math.  Anal., 
33  (2002),  pp.  1107  1116. 

[47]  T.  Lyche,  K.  Morken.  and  E.  Quak,  Theory  and  algorithms  for  non-uniform  spline 
wavelets ,  in  “Multivariate  Approximation  and  Applications”  (N.  Dyn,  D.  Leviatan,  D. 
Levin,  and  A.  Pinkus,  Eds.),  Cambridge  University  Press,  2001,  pp.  152  187. 

[48]  J.  Morgan  and  R.  Scott,  A  nodal  basis  for  C 1  piecewise  polynomials  of  degree  n  >  5, 
Math.  Comp.,  29  (1975),  pp.  736  740. 

[49]  G.  Niirnberger,  L.  L.  Schumaker,  and  F.  Zeilfelder,  Local  Lagrange  interpolation  by 
bivariate  C 1  cubic  splines ,  in  “Mathematical  Methods  in  CAGD:  Oslo  2000”  (T.  Lyche 
and  L.  L.  Schumaker,  Eds.),  Vanderbilt  University  Press,  2001,  pp.  393  404. 

[50]  P.  Oswald,  Hierarchical  conforming  finite  element  methods  for  the  biharmonic  equation , 
SIAM  J.  Nurner.  Anal.,  29  (1992),  pp.  1610  1625. 

[51]  P.  Oswald,  Multilevel  Finite  Element  Approximation ,  Teubner,  Stuttgart,  1994. 

[52]  J.  Peetre  and  G.  Sparr,  Interpolation  of  normed  Abelian  groups ,  Ann.  Mat.  Pura  Appl., 
92  (1972),  pp.  217  262. 

[53]  P.  Petrushev,  Direct  and  converse  theorems  for  rational  and  spline  approximation  and 
Besov  spaces ,  in  “Function  spaces  and  applications  (Lund,  1986)”,  Lecture  Notes  in 
Math.,  1302,  Springer,  Berlin,  1988,  pp.  363  377. 

[54]  P.  Petrushev,  Multivariate  n-term  rational  and  piecewise  polynomial  approximation ,  J. 
Approx.  Theory,  121  (2003),  pp.  158  197. 

[55]  P.  Petrushev  and  V.  Popov,  Rational  approximation  of  real  functions ,  Cambridge  Uni¬ 
versity  Press,  1987. 

[56]  P.  Sablonniere,  Quasi-interpolants  associated  with  H-splines  on  a  three  direction  mesh , 
J.  Comput.  Appl.  Math.,  66  (1996),  pp.  433  442. 

[57]  L.  L.  Schumaker,  On  supersplines  and  finite  elements ,  SIAM  J.  Numer.  Anal.,  26  (1989), 
pp.  997  1005. 

[58]  R.  Stevenson,  Piecewise  linear  (pre-)wavelets  on  non-uniform  meshes ,  in  “Multigrid 
Methods  V”  (W.  Hackbusch  and  G.  Witturn,  Eds.),  Springer,  Berlin,  1998,  pp.  306 
319. 


55 


Oleg  Davydov 
Mathematisches  Institut 
Universitat  Giessen 
D-35392  Giessen,  Germany 
Phone:  0641-99-32194 
Fax:  0641-99-32199 

E-mail:oleg .  davydovOmath .  uni-giessen .  de 

http://www.uni-giessen.de/www-Numerische-Mathematik/davydov/ 

Pencho  Petrusliev 
Department  of  Mathematics 
University  of  South  Carolina 
Columbia,  SC  29208 
Phone:  (803)  777-6686 
Fax:  (803)  777-6527 
E-mail:  pencho@math .  sc  .  edu 
http:  / / www. math. sc.edu/~pencho 


56 


